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PREFACE 



This text was prepared originally and has been revised annually by 
the Space Directorate of the Air University Institute of Professional 
Development. The Space Handbook serves as the text for the Funda- 
mentals of Space Operations Course, a resident course within AUIPD. 
As such, the text was written at an intermediate level of academic dif- 
ficulty but. with considerable depth of detail. 

The objectives of the Fundamentals of Space Operations Course are: 

1 . To provide the student with an understanding of the basic physical 
laws and principles of the space environment, propulsion, orbital 
mechanics, guidance and control, and atmosnherU penetration which 
permit and limit space operations. 

2. To provide the student familiarization with: objectives of the 
national space effort; current technology; propulsion devices and launch 
vehicles; electronic applications; present and possible future weapon and 
support systems including limitations and feasibility. 

3. To stimulate thought on new ideas and concepts so that the stu- 
dent may apply more effectively his knowledge in performance of space 
planning and operational duties. 

The Space Handbook also serves as the support text for the Astro- 
nautics and Space Operations phase of instruction, Air Force ROTC 
Aerospace Studies 300 Course. 

Recommendations for improvements of Space Handbook should be 
sent to: Commandant, Air University Institute for Professional De- 
velopment, Maxwell AFB, Alabama, 36112. 
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CHAPTER 1 



THE SPACE ENVIRONMENT 



F OR THOUSANDS of years man has looked at the heavens and wondered. 

What are the stars? What are the planets? The Moon? The shooting stars? 
Why and how do these heavenly bodies move? Answers derived from supersti- 
tion, philosophy, religion and fear abound in the literature and folklore of all 
peoples. Only recently, in the history of man, have answers been found in 
observation and experimentation. And. even these answers are tentative. The 
success of the first manned lunar landing, Apollo 11, was truly a milestone in 
the search for understanding which continues at a quickening pace. 

Today there are at least two ways of looking at the space environment. The 
first is the magnificent look— the look that sees space as the whole universe in 
terms of both matter and energy. The second is the practical look the look that 
sees space as another region in which man has begun travel. The former staggers 
the imagination and stimulates both wonder and reverence. The latter is the 
immediate concern of the military man. 

This chapter will start with a brief discussion of the universe as it is believed to 
be today. Following that will be a more detailed presentation of the characteristics 
of the near-earth space which has immediate military importance. 



THE UNIVERSE 

Man lives on an Earth which is one of nine planets in separate orbits around a 
star called the Sun. In addition to the planets, the Solar System contains thirty-two 
natural satellites, a variable number of artificial satellites, about thirty thousand 
asteroids, a hundred billion comets and countless specks of dust. These numbers 
seem impressive. But the Sun, which is the master of the entire system, is more 
impressive. It contains 99.9 percent of the matter in the Solar System. 

The nearest stellar neighbor to the Sun is a star called Alpha Centauri, a 

conspicuous double star which is visible from the southern hemisphere. About 
two degrees from it is Proxima Centauri, so named because it was once thought to 
be even nearer than Alpha Centauri. This group of three stars is about 4.3 light 

years from the Sun. That is, it takes light, traveling at 186,300 miles per second, 

about 4.3 years to reach the Solar System. It would take over 100,000 years for a 
spacecraft tiavelling at 25,000 miles per hour to make the trip. To date, the 
fastest man has ever travelled is about 25,000 miles per hour. It is clear that man 
cannot live long enough to travel to the nearest star. Of course, this is the present 
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status of such travel. In the future greater speeds, or perhaps the contraction of 
time which Dr. Albert Einstein predicted, may make such travel a reality. 

Another way to visualize these immense distances is to imagine the Sun as 
represented by a golf ball. On this scale Alpha Centauri would be represented by 
a pair ot golf balls which are a traction of a mile apart and 500 miles away. 
Proxima would be a grain of sand about 25 miles from the pair of golf balls. 
Some of the other closest stars are: 



Barnard’s Star 
Wolf 359 
Sirius 
Procyon 
Vega 



6.1 light years 
8 light years 

8.6 light years 
1 1 light years 
27 light, years 



But these are only the Sun’s nearest neighbors. Betelgeuse is 300 light years 
away. Polaris (the North Star) is 600 light years away. Yet all of these are only 
a few of a vast array of stars that form a group called the Milky Way. 

The name given to a large group of stars, dust, and gas that stay together in a 
structure is a galaxy. The Milky Way is simply the view from Earth of the 
galaxy in which the Sun is one star. Figure 1 is a view of the Milky Way as it 
would appear from the side. Figure 2 is its appearance from above. Of course 
these are artists’ conceptions, because such pictures cannot be taken from inside 
the galaxy. The Milky Way is about 100,000 light years in diameter. The Solar 
System is located in one of the spiral arms of the galaxy, about 30,000 light years 
from the center. 

How many stars are there in this magnificent structure? Of course no one 
knows the e.x .ct answer. However, Dr. Harlow Shapley estimates the number at 
about 200 .ion. All of these arc in motion around the center of the galaxy. At 

the distar, of the Solar System from the center of the galaxy, the speed of 
revolution is about 135 miles per second or about 486,000 miles per hour. Even 
at this tremendous speed, it takes the sun 220 million years to make one trip 
around .he galactic center. Stars closer to the center move faster and those 
farther from the center move more slowly. 

Beyond the Milky Way the telescopes show other objects. What are they? Gas? 
Dust? Stars? In the year 1755 Immanuel Kant suggested that these were “island 
universes” — other galaxies similar to the Milky Way, each consisting of billions of 
stars. However, it was not until 1917 that an astronomer using the Mount Wilson 
telescope identified a star in one of these objects beyond the Milky Way. 

Today we know that these objects are indeed other galaxies. Seme have 
diameters in the order of 7,000 light years. Others have diameters even greater 
than the Milky Way — about 150,000 light years. The smaller ones probably 
contain a billion stars. The larger ones may contain two hundred billion or more. 
How many galaxies are there in the universe? Again the answer must be an 
estimate, and the most recent one is huge — one hundred billion! Finally, if one 
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Figure 1 
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nsks the question, “How many stars are there in the Universe?” the answer is 
approximately 10- 1 stars! Although no one can imagine this number, Sir James 
Jeans has provided an analogy which helps. He suggests that the number of stars 
in the Universe is something like all the grains of sand on all the beaches of all 
the Earth. 

To complete the modern idea of this stupendous structure, it is believed that the 
whole thing is expanding. Every galaxy is rushing away from every other galaxy 
at tremendous speed! The most distant galaxy visible in a telescope is racing 
away at about 75,000 miles per second or 270 million miles per hour. 

THE MILITARY SPACE ENVIRONMENT 

Most of the space discussed above is of little military importance today. How- 
ever, the near-earth space is important for many military operations. Its charac- 
teristics set boundary conditions on both the operations that can be conducted 
and the equipment required. Consequently, the remainder of this chapter will be 
devoted to the hazards and physical conditions of near-earth space. “Near-earth” 
will not be defined precisely, but it is understood to mean not more than 100 
million miles from Earth. The Sun is about 93 million miles away. Thus the Sun, 
the inner planets, and the space between them is the concern of this section. 
Emphasis will be upon space this side of the Moon. 

The first topic to consider is: Where does space begin? Certainly it does not 
begin at the surface of the Earth because that is where the atmosphere begins. At 
an altitude of 10,000 feet the oxygen pressure of the atmosphere is not great 
enough to keep people efficient over a long period of time. Many people, of 
course, become acclimatized to altitudes of 10,000 feet and higher. But for a man 
who lives near sea-level, the oxygen pressure at levels above 10,000 feet is 
insufficient to sustain active and efficient performance. Thus the Air Force 
requires the use of supplemental oxygen by crew members at altitudes above 
10,000 feet. 

Approximately one half of the Earth’s atmosphere is below an altitude of 
three miles. But it is not until an altitude of about nine miles that supplemental 
oxygen fails as a sufficient aid to sustain human life. Here the combined pressure 
of carbon dioxide and water vapor in the lungs equals the outside atmosphere 
pressure and breathing cannot take place without supplemental pressure. Hence, 
at this altitude pressure cabins or pressure suits become a necessity. 

The vapor pressure of man’s body fluids is about 47 mm of mercury. As soon 
as the atmospheric pressure drops to this level, bubbles of water vapor and other 
gases appear in the body fluids. This means literally that the blood will boil. 
The gas bubbles first appear on the mucous membranes of the mouth and eyes 
and later in the veins and arteries. This would happen to an unprotected man at 
an altitude of 12 miles. However, supplemental pressure will suppress this evil. 

At 15 miles compressing the outside air to pressurize a cabin is no longer 
effective. At that altitude the air density is about % 7 of the sea level value. 
Compressing this thin air is perhaps not an impossible task but it certainly is an 
uneconomical task. Further, the act of compressing air would involve undesirable 
heat transfer to the air. Finally, at this altitude the atmosphere contains a 



significant percentage of ozone. If this were compressed it would poison the cabin 
atmosphere. Hence, above this altitude the cabin or space suit must have a 
supply of both oxygen and pressure independent of the outside atmosphere. As 
far as man is concerned, 15 miles above the surface can be considered the 
beginning of space. Above this altitude man must take everything he needs with 
him. His environment will supply him with neither food nor air. He needs a 
seated environment containing necessary supplies from Earth. 

Five miles farther out, at the 20 mile level, is the operating limit for turbojet 
engines. At 28 miles ramjet engines do not have enough air to operate. Above 
this altitude engines must be supplied with both a fuel and an oxidizer. Thus, to a 
propulsion engineer 28 miles above the Earth is the beginning of space. Above 
this he must use rockets. 

In one sense space begins at 50 miles because flight above this altitude earns a 
crew member the right to wear Astronaut’s wings. 

In 1964, a New York law firm asked the Air Force Office of Aerospace 
Research to define the beginning of space. This scientific organization based the 
answer on aerodynamic forces. Such forces acting on ballistic reentry vehicles, 
lifting reentry vehicles and boost-glide orbital vehicles can usually be neglected 
at altitudes above 100 kilometers or 62 miles. Thus, for the aeronautical engineer 
concerned with lift and drag, space may begin at 62 miles. 

At about 100 miles above the earth is a region of darkness and utter silence. 
This is the region of the black sky. The stars appear as brilliant points of 
light, and between them is absolute black because there is not enough air to 
scatter light. Neither is there enough air to carry sound or shock waves. There 

are no sonic booms. 

From the above discussion it is clear that there are many answers to the 
question, “Where does space begin?” The acceptable answer depends upon the 
reference frame in which the question was asked. 

“Is space really nothing?” The answer is “No.” Space is filled with surpris- 
ing amounts of matter and is flooded with energy. First consider the density of 
matter by starting up from the Earth’s surface. At the surface the concentra- 
tion of particles in air is about a million, million, million particles per cubic 
centimeter (10 18 /cm 3 ). There is a decrease in particle density with altitudes ard 
the average figures given are only approximate. An average figure for the zoik-j 
between 7 miles and 50 miles is about 10 14 /cm 3 . From 50 miles to 600 miles 
an average figure is about a million particles per cubic centimeter (10 6 /cm 8 ). 
From 600 miles to 1200 miles there ere still in the order of 100 particles per 

cubic centimeter. Above 1200 miles there will be found something like one 

particle per cubic centimeter. This is certainly far from nothing. There are also 
localized conditions that cause the particle density to be much higher. Some o 
these conditions will be discussed later. And, all of space is flooded with electro- 
magnetic energy in many forms from the Sun, from the stars in the Milky 
Way and even from other galaxies in the Universe. 

However, from another point of view, space is “not much. Consider air 

pressure. At the Earth’s surface the pressure varies around the figure 760 mm 

of mercury. Above 1200 miles the pressure is much less than one mm of mer- 
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cury. In fact, it probably is around 10~ 12 to 10~ 1C Torr.* This pressure is so 
low that it is often called a “hard vacuum.” It causes some unexpected phe- 
nomena. A few of these will be discussed as illustrations. 

In the atmosphere any metal is covered with at least a single layer of ab- 
sorbed gas. In a hard vacuum this film of gas bleeds into space. Metals touch- 
ing each other tend to weld together. In the atmosphere this doesn’t happen 
because the thin film of air acts as a lubricant keeping the metals apart. To 
prevent this “cold welding” in space, special measures must be taken. 

Some metals are stronger in a hard vacuum. If a crack forms in a metallic 
surface when the metal is surrounded by air, molecules of air immediately enter 
the crack. Chemical reaction with the metal occurs. If the reaction product is 
more voluminous than the original crack, a wedging action occurs and enlarges 
the crack. In a hard vacuum a chemical reaction causing an enlargement of a 
crack does not occur. Thus some metals may be stronger in space than they 
are on Earth. 

To study space effects such as the above, it is useful to simulate a hard 
vacuum on Earth. Late in 1965 USAF completed a large Aerospace Environ- 
mental Simulation Chamber at the Arnold Engineering Development Center in 
Tennessee. 



Electromagnetic Radiation 

Visible light, ultra-violet, x-rays, infrared, radio and other forms of energy 
can travel through the hard vacuum of space as electromagnetic radiation. This 
term refers to the fact that radiation consists of a varying electric field and a 
varying magnetic field. Together these fields form a wave. Such a wave can 
be transmitted through a vacuum and does not require the presence of a ma- 

* A T°rr is the same as one millimeter of mercury. It is named in honor of Torricelli and is commonly 
used in low pressure measurements. 




Figure 3 



terial Medium, This fo m of energy floods all of space. Its intensity varies 
with proximity to the Sun, or to a star. 

Notice in Figure 3 that visible light covers only a very narrow band in the 
spectrum of electromagnetic energy, The entire spectrum ranges from frequen- 
cies of about 10* Hertz up to about 10-' Hertz. 

Most of the radiation coming from the Sun and stars is absorbed by the 
Earth’s atmosphere before it reaches the surface of the Earth. In fact, there 
are only two “windows” through which space may be observed. One window 
includes visible light frequencies and part of the ultraviolet and infrared fre- 
quencies. This is called the optical window. Another window called the radio 
window is found in the radio frequencies of approximately 10 9 cycles per sec- 
ond. Man is protected from (or denied, depending upon the point of view) all 
other frequencies by the Earth s atmosphere. In space radiation of frequencies 
throughout the spectrum is present. Man may use it as a source of informa- 
tion, but he must also protect himself against it by a space cabin or space suit. 

Meteoroids and Micrometeoroids 

The three terms “meteoroid,” “meteor,” and “meteorite” have similar meaning 
and are often used interchangeably. Meteoroid refers to a particle, large or 
small, moving in space. When a meteoroid enters the atmosphere and begins 
to glow, it is then called a meteor. If that same particle survives the trip through 
the atmosphere and hits the earth, the remnant is called a meteorite. Some 
meteoroids must be very large because meteorites with masses of several tons 
have been found. Most meteoroids, however, are quite small. Extremely small 
meteoroids are called micrometcoroids. 

Meteoroids and micrometcoroids move with speeds varying from about 30,000 
miles per hour to 160,000 miles per hour. At these speeds, impact between a 
satellite and a large meteoroid would be catastrophic. Impacts between micro- 
meteoroids and a satellite would not be catastrophic but could erode the satel- 
lite’s surface. 

Several satellites have rather mysteriously ceased to function, and there is 
some conjecture that meteoroid impact damage may have been the cause. The 
probability of meteoroid and micromctoroid impact has been extensively studied. 
Many methods have been employed. Micrometeoroids so small that it would 
take about 125 of them to equal the thickness of a piece of paper have been 
captured, and they have been studied with an electron microscope. The Pegasus 
satellites, Explorer XXIII and Explorer XVI, are examples of satellites used to 
study the problem. The Pegasus satellites reported the number of penetrations 
of their panel materials such as aluminum. Explorer XVI, launched December 
16, 1962, reported 62 meteoroid penetrations in I-V 2 months of space travel. 
Mariner IV was hit by over 1 50 rnicrometeoroids on its trip to Mars. These 
and other data have led to the conclusion that the probability of a satellite 
being hit by very small micrometcoroids is high and of being hit by catastrophi- 
cally large ones is small. 

Figure 4 provides some idea of the probabilities involved in the Apollo 
Project. The probabilities are based upon an estimated vehicle cross section of 
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Figure 4 

ten square meters and an exposure of ten days. The horizontal axis shows mass 
of the meteoroids or micrometeoroids in grams. (About 5 grams equals the 
weight of a nickel.) The vertical axis is scaled in probability until the figure 
one is reached. This figure means that the vehicle is certain to be hit. Above 
that, the scale means nun 1 her of hits in a ten-day period. Notice that the 
probability ol collision with a particle which weighs 10^ r> grams is 1.0. This 
kind of particle would penetrate about 0.05 cm into an aluminum skin. 

In general, the problem of meteoroid hazard is not as great as was once 
thought, but it is not negligible. Protection is provided in space suits and cap- 
sules. The possibility of a catastrophic hit remains. Extensive studies are con- 
tinuing. 



Cosmic Rays 

Cosmic rays are very small particles which have been travelling within the 
Milky Way for millions of years. Some may even come from other galaxies. 
They arc positively charged and move with great speed. It would take about 10 
billion cosmic rays to equal the thickness of a piece of paper. About 84% of 
them are protons, and about 14% are Alpha particles. The rest are nucleii of 
atoms with higher atomic numbers ranging from lithium to iron. The speed 
of cosmic rays is not usually defined. Rather, the kinetic energy level is given. 
Speed and mass are the two factors which are incorporated into kinetic energy. 
The energy of cosmic rays is usually stated in electron volts (ev).* In these 

* Other frequently used terms are kev (one thousand ev). mev (one million ev) and bev (one billion ev). 
Energy in electron volts can be converted to kinetic energy. Then if the mass of a proton is known, its 
speed can be calculated. 



units, the encrcv of cosmic rays varies from about 40 million electron volts to 
about a million' million, million electron volts (10>- million ev). Energies greater 
than a billion electron volts mean that the particle's speed is neai 1 that of light. 

Cosmic rays arrive at Earth from all directions. Additional piotons also arrive 
at Earth from the direction of the Sun. 1 he former are otten called galactic 
cosmic rays"; the latter are "solar piotons." When cither galactic cosmic rays or 
solar protons collide with atoms of the Earth’s atmosphere, the atoms break up, 
or disintegrate. This process produces a shower of smaller particles which me ve 
off in all directions from the site of the collision. The particles are varied in 
kind and energy but include neutrons, electrons, and mesons, as well as some 
heavier particles. One of the ways used to study cosmic radiation and its 'ollision 
with atoms is to fly a photographic plate into space. When the plate is recovered, 
it frequently will show the result of a cosmic ray striking an atom in the photo- 
graphic plate and breaking the atom into smaller particles. The charged particles 
leave traces on the plate; the resulting image resembles a star. Such a collision 
is often called a “collision star'’ or “disintegration star.” 

The Earth’s atmosphere serves to protect man from the effects of primary 
cosmic radiation; in space man would be exposed to a higher intensity of cosmic 
radiation than he is on Earth. Also, since the cosmic radiation particles are 
charged particles the Earth’s magnetic field acts as a shield. The effect is to 
divert many of Hie particles to the north and to the south. The result is that 
man on earth is protected by both the Earth s atmosphere and magnetic field. 
When man ventures into space, he abandons this natural protection. 

Beyond the magnetic influence of the Earth, the cosmic radiation from the 
galaxy presents a flux of one to two particles per square centimeter per second. 
Despite the high energy of galactic cosmic rays, there will be little hazard to 
man because of the low flux. Even in a year’s time the total radiation dose 
expected from this source is from 6 to 20 rads. This is less than one-tenth of the 
dose that would make one-half of recipient humans feel sick or nauseated. 

Solar Flares 



The high speed solar protons emitted by a solar flare are probably the most 
potent of the radiation hazards to space flight. Flares themselves are probably the 
most spectacular disturbances seen on the sun. They are observed optically as 
a sudden, large increase in light from a portion of the Sun’s atmosphere. A flare 
may spread in area during its lifetime which may be from several minutes to a 
few hours. Flares are classified according to a range of importance of 0 to 4. 

There is a relationship between the number of sunspots and the frequency of 
flare formation, but the most important flares do not necessarily occur at sun- 
spot maximum. 

There are many events that may occur at Earth following a solar flare although 
not many flares produce all of the possible events. In addition to the increase in 
visible light, minutes after the start of a flare there is a Sudden Ionospheric Dis- 
turbance (STD) in the Earth’s ionosphere. This, in turn, causes short wave fade- 
out, resulting in the loss of long range communications for 15 minutes to 1 hour. 
X-rays emitted by the flare probably cause the SID’s. During the first few minutes 
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Table 1 

Solar Flare Classification 





Area* 


* — 


Importance 


(Millionths of Solar Hemisphere ) A verage Duration 


6 


Less than 100 


17 min 


i 


100 to 249 


32 min 


2 


250 to 599 


69 min 


3 


600 to 1200 


145 min 


4 


Greater than 1200 


145 min 


BRIGHTNESS CATEGORIES FAINT (F) 






NORMAL (N) 
BRILLIANT (B) 




* The area of the earth’s disk is approximately equivalent to the area of an Importance one flare 


of a flare there may be 


a radio noise storm, consisting of bursts 


.T noise over a 



wide range of frequencies. In addition, there may be disturbances in th- Earth’s 
magnetic field, changes in the Auroras, and decreases in galactic cosmic ray in- 
tensity. However, from the point of view of a space traveller, by far tr.e most 
important effect is the marked increase in solar protons. The energy of these 
protons ranges from about 10 million electron volts to about 500 million electron 
volts. The flux may be quite high. Consequently, the dose of radiation accumu- 
lated during exposure to the solar protons may vary from negligible to well above 
a lethal dose. 

Construction of space vehicles sufficiently shielded to protect against all possi- 
ble solar proton events is impractical due to the amount and density of required 
shielding materials. Consequently the best hope for protection lies in developing 
a method for reliable prediction of flare occurrence and intensity. The USAF Air 
Weather Service is charged with this responsibility. 

The Van Allen Radiation Belt 

Another problem that man must overcome in venturing into space is trapped 
radiation. As a result of experiments conducted in 1958 with the US Explorer 
I satellite and subsequent experiments, Dr. James A. Van Allen and his asso- 
ciates discovered the existence of geomagnetically trapped particles encircling the 
earth. When electrons and protons, and perhaps some other charged particles, 
encounter the Earth’s magnetic field, many of them are trapped by the field. 
They oscillate back and forth along the lines of force, and since the magnetic field 
completely encircles the Earth, the trapped particles completely encircle the Earth. 

The belt has an inner and outer portion. Recent data shows that the toroidal 
shaped volume occupied by the Van Allen radiation is permeated with both pro- 
tons and electrons. The protons are most intense at about 2,200 miles. The elec- 
tron flux peaks at about 9,900 miles. The low particle density separating the two 
belts is often called the “slot.” In this particular volume of space some phenomena, 
as yet not fully undestood, reduces the lifetime of the charged particles. 

The inner Van Allen belt starts at an altitude of about 250 miles to 750 miles, 
dependiilg upon latitude. It extends to about 6,200 miles where it begins to over- 
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lap the outer belt and where the “slot’’ begins. The inner belt extends from about 
45 north latitude to about 45' south latitude. 

The outer Van Allen belt begins at about 6,200 miles and extends to between 
37.000 and 52,000 miles. The upper boundary is dependent upon the activity of 
the Sun. 

Both the inner and outer belts were affected by a high altitude nuclear device 
test in July 1962. Radiation in both belts increased after the detonation, and the 
low radiation slot separating the two belts was eliminated for some time. 

Experience has shown that space vehicles in low circular orbit (125-350 miles) 
receive an insignificant amount of radiation from the Van Allen zones. However, 
a vehicle in a highly eccentric orbit or one in a high altitude circular orbit can re- 
ceive an important dose. For example, a satellite in a synchronous orbit over the 
equator will be close enough to the center of the outer zone to accumulate a 
hazardous dose. But, as demonstrated by Apollo lunar missions, man can safely 
transit these zones in a spacecraft with minimal shielding by judicious selection of 
the flight trajectory. 

The Van Allen belt varies daily with changes in the magnetosphere. On the 
Sun side of Earth it is flattened. On the night side of the Earth it is elongated. 

The Solar Wind 

Because of the high temperature of the Sun’s corona, protons and electrons 
beyond a certain distance from the San acquire velocities in excess of the escape 
velocity from the Sun. Thus there is a continuous outward flow of charged par- 
ticles in all directions from the Sun. This has v ,cn called the solar wind. It is a 
plasma wind, rather than a gas wind. Its velocity and density vary with sunspot 
activity. During the time of sunspot minimum at the Earth’s distance from the 
Sun, the density is about 100 particles per cubic centimeter. The speed is about 
300 miles per second. At sunspot maximum the corresponding density is probably 
around 10,000 particles per cubic centimeter, and the speed is about 900 miles 
per second. When the solar wind encounters the Earth’s magnetosphere, it flows 
around the magnetosphere, which is flattened in the process. On the dark side of 
the Earth the magnetosphere is elongated. 

The energy of the particles in the ;:olar wind is not high. No hazard to man 
is expected from the wind. However, it is either the cause of or a contributor to 
the aurora which illuminates the polar night sky. Within the past few years ground 
based observations h.s s been combined with information acquired by rockets 
and earth satellites ‘o provide an explanation of this previously baffling beauty. 
The magnetosphere of the Earth acts like a giant cathode-ray tube, directing 
charged particles from the solar wind into beams and focusing them on the Earth’s 
polar regions. The aurora is a fluorescent luminosity produced by the electrons 
and protons of the solar wind which strike atoms and molecules of oxygen and 
nitrogen high in the atmosphere of ti e polar region. 

The radiation hazard of space may be summarized as follows: 

( 1 ) Electromagnetic energy in space may be shielded and is not a serious threat 
to life. 
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(2) Electromagnetic energy may upset radio communication and guidance 
equipment. 

(3) Galactic cosmic rays and solar wind do not present a serious threat to 
space travel. 

(4) Van Allen radiation does present a serious threat, but the location of the 
belts is sufficiently well known that flight trajectories may be planned to 
limit time spent in the hazardous regions. 

(5) Protons emitted at the time of a solar flare present the greatest uncertainty 
and the greatest threat to manned flight in regions beyond the protection 
of the Earth’s atmosphere and magnetosphere. 
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CHAPTER 2 



ORBITAL MECHANICS 



T HE STUDY of trajectories and orbits of vehicles in space is not a new science 
but is the application of the concepts of celestial mechanics to space vehicles. 
Celestial mechanics, which is mainly concerned with the determination of tra- 
jectories and orbits in space, has been of interest to man for a long time. When 
the orbiting bodies are man-made (rather than celestial), the topic is generally 
known as orbital mechanics. 

The early Greeks postulated a fixed earth with the planets and other celestial 
bodies moving around the earth, a geocentric universe. About 300 B. C., Aristar- 
chus of Samos suggested that the sun was fixed and that the planets, including 
the earth, were in circular orbits around the sun. Because Aristarchus’ ideas were 
too revolutionary for his day and age, they were rejected, and the geocentric 
theory continued to be the accepted theory. In the second century A.D., Ptolemy 
amplified the geocentric theory by explaining the apparent motion of the planets 
by a “wheel inside a wheel” arrangement. According to this theory, the planets 
revolve about imaginary planets, which in turn revolve around the earth. It is 
surprising to note that, even though Ptolemy considered the system as geocentric, 
his calculations of the distance to the moon were in error by only 2%. Finally, in 
the year 1543, some 1800 years after Aristarchus had proposed a heliocentric 
(sun-centered) system, a Polish monk named Copernicus published his De 
Revolutionibus Orbium Coelestium, which again proposed the heliocentric theory. 
This work represented an advance, but there were still some inaccuracies in the 
theory. For example, Copernicus thought that the orbital paths of all planets were 
circles and that the centers of the circles were displaced from the center of the 
sun. 

The next step in the field of celestial mechanics was a giant one made by a 
German astronomer, Johannes Kepler (1571-1630). After analyzing the data 
from his own observations and those of the Danish astronomer Tycho Brahe, 
Kepler stated his three laws of planetary motion. 

A contemporary of Kepler’s, named Galileo, proposed some new ideas and 
conducted experiments, the results of which finally caused acceptance of the 
heliocentric theory. Some of Galileo’s ideas were expanded and improved by 
Newton and became the foundation for Newton’s three laws of motion. Newton’s 
laws of motion, with his law of universal gravitation, made it possible to prove 
mathematically that Kepler’s laws of planetary motion are valid. 

Kepler’s and Newton’s work brought celestial mechanics to its modern state 
of development, and the major improvements since the days of Newton have 
been mainly in mathematical techniques, which make orbital calculations easier. 
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Because the computation of orbits and trajectories is the basis for predicting 
and controlling the motion of all bodies in space, this chapter describes the fun- 
damental principles of orbital mechanics upon which these computations are 
based. It also shows how these principles apply to the orbits and trajectories used 
in space operations. 



MOTION OF BODIES IN ORBIT 

Bodies in space move in accordance with defined physical laws. Analysis of 
orbital paths is accomplished by applying these laws to specific cases. Orbital 
motion is different from motion on the surface of the earth; however, many con- 
cepts and terms arc transferable, and similar logic can be applied in both cases. 
An understanding of simplified linear and angular motion will permit a more 
thorough appreciation of a satellite’s path in space. 

Linear Motion 

Bodies in space are observed to be continuously in motion because they are 
in different positions at different times. In describing motion, it is important to 
use a reference system. Otherwise, misunderstanding and inaccuracies are likely 
to result. For example, a passenger on an airliner may say that the stewardess 
moves up the aisle at a rate of about 5 ft per sec, but, to the man on the ground, 
the stewardess moves at a rate of the aircraft’s velocity plus 5 ft per sec. The 
man in the air and the man on the ground are not using the same reference 
system. For the present, the matter of a reference system will be simplified by 
first describing movement along a straight line, or what is called rectilinear 
motion. 

Rectilinear motion can be described in terms of speed, time, and distance. 
Speed is the distance traveled in a unit of time, or the time rate of change of 
distance. An object has uniform speed when it moves over equal distances in 
equal periods of time. Speed does not, however, completely describe motion. 

Motion is more adequately described if a direction as well as a speed is given. 
A speedometer tells how fast an automobile is going. If a direction is associated 
with speed, the motion is now described as a velocity. A velocity has both a 
magnitude (speed) and a direction, and it is therefore a vector quantity. 

Uniform speed in a straight line is not the same as uniform speed along a 
curve. If a body has uniform motion along a straight line for a given time, then 

average velocity is represented by the equation v = ~ — ;° . In the equation, 

to 

s 0 is the initial position, s* is the final position, t 0 is the initial time, and t f is the 
final time; or more simply, the velocity is the change in position divided by the 
change in time. The units of velocity are distance divided by time, such as ft 
per sec or knots (nautical miles per hour).* Since velocity is a vector quantity, 
it may be treated mathematically or graphically as a vector. 

If velocity is not constant from point to point (i.e., if either direction or speed 
is changed), there is acceleration. Acceleration, which is also a vector quantity, 

• The nautical mile (NM) is one minute of a great circle. In this course, use the conversion that 1 NM = 
6,080 feet =1.15 statute miles. 



is the time rate of change of velocity. The simplest type of acceleration is one 
in which the motion is always in the same direction and the velocity changes 
equal amounts in equal lengths of time. If this occurs, the acceleration is constant, 
and the motion can be described as being uniformly accelerated. 

The equation a nv = ■ _ I 0 - defines the average acceleration, over the 

tf to 

specified time interval. A good example of a constant acceleration is that of a 
free-falling body in a vacuum near the surface of the earth. This acceleration 
has been measured as approximately 32.2 ft per sec per sec, or 32.2 ft per sec J . 
It is usually given the symbol g. Since an acceleration is a change in velocity over 
a period of time, its units are ft/sec 2 , or more generally, a length over a time 
squared. Actually, constant acceleration rarely exists, but the concepts of constant 
acceleration can be adapted to situations where the acceleration is not constant. 

The following three equations are useful in the solutions of problems involving 
linear motion: 



where s is linear displacement, v„ is initial linear velocity, v f is final linear velocity, 
a is constant linear acceleration, and t is the time interval. 



If a particle moves along the circumference of a circle with a constant tangential 
speed, the particle is in uniform circular motion. Since velocity signifies both 
speed and direction, however, the velocity is constantly changing because the di- 
rection of motion is constantly changing. Now, acceleration is defined as the time 
rate of change of velocity. Since the velocity in uniform circular motion is chang- 
ing, there must be an acceleration. If this acceleration acted in the direction of 
motion, that is, the tangential direction, the magnitude of the velocity (the speed) 
would change. But, since the original statement assumed that the speed was 
constant, the acceleration in the tangent'al direction must be equal to zero. There- 
fore, any acceleration that exists must be perpendicular to the tangential direction, 
or in other words, any acceleration must be in the radial direction (along the 
radius). 

Average speed is equal to the distance traveled divided by the elapsed time. 
For uniform circular motion, the distance in one lap around the circle is 2 ttt , 
which is covered in one period (P). Period is the time reqmred to make one 
trip around the circumference of the circle. Therefore, the tangential speed 

v t = ^ b F . In uniform circular motion, the particle stays the same distance from 

the center, therefore, radial speed, v r = 0. It has already been shown that 
a t = 0; and it will be shown in the next section that a r = 



(1) s 

(2) v f 

(3) 2as 




Angular Motion 



r 



Translatory motion is concerned with linear displacement, s; velocity, v; and ac- 
celeration, a. Angular motion uses an analogous set of quantities called angular dis- 
placement, 6; angular velocity, co ; and angular acceleration, a. 

In describing angular motion, it is convenient *o think of it in terms of the rota- 
tion of a radius arm (r), as shown in Fi^uie 1. The radius arm initially coincided 
with the polar axis, but at some time later (t seconds) it was positioned as shown. 




Figure 1. Position of radius arm as rotated one radian (57.3 degrees) from the starting 

point. 

Angular displacement (9) is measured in degrees or radians. A radian is the 
angle at the center of a circular arc which subtends an arc length equal to the radius 
length. If the length of s equaled the length of r, 9 would be equal to one radian, or 
57.3°. The central angle of a complete circle is 360° or 2i r radians (27 r = 6.28). 

The following equations for angular motion are analogous to those studied earlier 
for rectilinear motion: 



6 — —^—radians ^ 

r _ * 

a> av = ^-rad/sec 

t 0 

«av = ^-rad/sec 2 

if to 

cot = co 0 + at 
a - - at 2 

6 — 0) o + — 2 — 

2 ad = o) t 2 — o) 0 2 

In the equations, 6 t is final angular position; 9„ is initial angular position; s is linear 
displacement (arc length); r is the radius; co is average angular speed; o) t is final 
angular speed; a> 0 is initial angular speed; t f is final time; t 0 is initial time; and a is 
constant angular acceleration. 

If a body is rotating about a center on a radius r, the tangential linear quantities 
are related to the angular quantities by the following formulas [where 6 , a), and a 
are in radians] : 

s = r$ 
v t = TCO 
a t = roc 
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Principles of the Calculus Applied to Astronautics 

Computations in the calculus are based upon the idea of a limit of a variable. Ac- 
cording to the formal definition, the variable x is said to approach the constant 1 as 
a limit when the successive values of x are such that the absolute value of the differ- 
ence x — 1 ultimately becomes and remains less than any preassigned positive num- 
ber, however small. 

An example will make the definition easier to understand. The area of a regular 
polygon inscribed in a circle approaches the area of the circle as a limit as the num- 
ber of sides of the polygon approaches infinity (Fig. 2). 





Figure 2. Increase in the number of sides of a regular polygon inscribed in a circle. 

The area of a triangle is 1/2 bh. In general, if there are n sides to a polygon, the 
polygon is made up of n triangles as shown in Figure 2. Therefore, the area of the 
polygon is 1/2 nbh. As the number of sides (n) approaches infinity as a limit, the 
product nb approaches the circumference of the circle (c). Also, as n approaches in- 
finity, the value of h approaches the radius (r) as a limit. 

lim 1/2 nbh = cr 
n oo 2 

This is read, “The limit of 1/2 nbh as n approaches infinity is equal toH.” 

But c = 2-7T r 

lim area of the polygon = lim nbh _ (27rr)r 
n co n-»oo2 2 

and ( 27rr )L = «7jT 2 = area of the circle 
2 

Now, an increment is the difference in two values of a variable. In the example 
above, the ’".crease in area when the inscribed polygon increases the number of sides 
by one is an increment of area; that is, the area of an inscribed square minus the 
area of an inscribed triangle is an increment of area. An increment is written as Ax 
which is read “delta x,” and does not mean A multiplied by x. 
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In the previous section the radial acceleration for uniform circular motion was 
given as a r = — . With the concepts of an increment and a limit, the value for 

radial acceleration can be determined mathematically. In Figure 3. an increment of 
arc has been expanded to permit closer examination. The length r is the distance 
from the center of the circle to the circumference. The horizontal distance v,At is 




Figure 3. An increment of an arc (left) and the increment expanded (right) to show 

change in velocity. 



the distance a body in uniform motion with a velocity v t would move in the time At. 
However, at the completion of the increment of time the body is not at point A but 
at point B, because this is uniform circular motion. The distance from A to B is 

equal to v r A t + ar - where the subscript r refers to radial. However, for uni- 



form circular motion v r = 0. Therefore, the distance AB is equal to 
applying the Pythagorean theorem to the triangle, 



a r A t- 
2 



. Now, 



r 2 + (v, At) 2 = £ r + -iftLj 

or, r 2 + v t 2 A t 2 = r 2 + r a r A t 2 + Ur - - — 

4 



Subtracting r 2 from both sides, 



v t 2 At 2 = r a r A t 2 + — 



Dividing both sides by A t 2 



v t 2 



r a r + 



a r 2 A t 2 

4 
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To find the instantaneous vaiues take the limit as A t -> 0. 



v, 2 = r a, + ar8 4 (0) = r a, 

I) 

a r = As was to be demonstrated. 



This text does not attempt to teach the processes of differentiating and integrating, 
but its purpose is to give the student some understanding of how the calculus is used 
in the study of space. 

The definition of the derivative of y with respect to x is, in symbol form, 

_ lim _ A _ y _ Any ca i cu i us book has a table of derivatives, and (here is also 
Ax-»0 A x 

one in The Engineer’s Manual by Ralph G. Hudson on pages 31 and 32. 

The average velocity over a period of time, as given in the previous section, is: 



v = Sf s ° 
tft 0 

Usually the average velocity is not of direct value in analysis, but the instanta- 
neous velocity is. The speedometer in a car measures instantaneous speed, and if a 
motorist is arrested for speeding, it is because of his instantaneous velocity, not his 
average velocity. If s is the path of a particle, its instantaneous velocity is equal to: 
ds = lim A s 

dt At -» 0 A t ‘ 

Example : A particle moves so that its distance from the origin at any time follows 
the formula s = t 3 . Find its average and final, velocity and acceleration after 3 sec- 
onds. 



y ss Sf S 0 to 0, S 0 0 

^ _t o tf =3, Sf = 27 

V a y = ^ 2—0 = ^ Answer 

V = _ds_ = _d_ ( t3 ) 

f dt dt 



From page 32 of The Engineer’s Manual: 



_d_ 

dx 



(u n ) 



n u n_1 



du 

dx 



d f+8'4 — 3t 8 "i dt = ^t 2 

dt 1 ; dt 

Vf = 3(3) 2 = 27 Answer 

a av = y° = 21 - f ° - = 9 Answer 
tf t 0 3 

a f = = 2 (3 t2_1 )“^ _ = 6t = IS Answer 
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Note that with the use of differential calculus, final or instantaneous values for ve- 
locity and acceleration can be determined, but only average values can be deter- 
mined from the formulas given in the previous section. 

If, in the example above, the acceleration were given as a f = 6t. the instanta- 
neous velocity and position could be determined by the process of integration. Inte- 
gral calculus is a summation process that is the inverse of differential calculus. 

Example: a f — 6t. Find v, after 3 seconds. If a curve is drawn with acceleration 
on the vertical axis and time on the horizontal axis, the area under the curve is the 
velocity (Fig. 4). Integration gives the sum of all the individual shaded rectangles as 



A t approaches 0 as a limit. As A t~*0, the area of the rectangles approaches the area 
under the curve as a limit and is the velocity in this problem. The symbol for inte- 
gration is f . From the table of fundamental theorems on integrals (Hudson’s The 
Engineer’s Manual, p. 39), 



In this example problem, the limits of integration, t — 0 to t — 3, are specified, so 
the + c (constant of integration) may be dropped. 




At 



Figure 4. Graph of the continuous function af — 6t. 





O 0 



Vf = 3t 2 I 3 = 3 (3) 2 - 3 (0) 2 = 27 Answer 



O 
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The processes of integration and differentiation of variables as applied to the 
computation of velocity and acceleration through the calculus are part of the study 
of motion taken up in the branch of dynamics known as kinematics. The study of 
the forces causing the motion belongs to another branch of dynamics called kinetics. 



LAWS OF MOTION 

Natural bodies in space follow the basic laws of dynamics, as described by 
Newton’s universal law of gravitation and his three laws of motion. By applying 
the basic laws and making use of calculus (also developed by Newton), one can 
explain and prove Kepler’s three laws of planetary motion. It would be well to 
review Kepler’s laws before stating Newton’s law of universal gravitation, which 
is one of the laws upoa which computation of trajectories and orbits* is based, 
and Newton’s three laws of motion, which describe terrestrial motion as well as 
celestial mechanics. 

Kepler’s Laws 

From his observations and study, Kepler concluded that the planets travel 
around the sun in an orbit that is not quite circular. He stated his first law thus: 
The orbit of each planet is an ellipse with the sun at one focus. 

Later Newton found that certain refinements had to be made to Kepler’s first 
law to take into* account perturbing influences. As the law is applied to manmade 
satellites, we must assume that perturbing influences like air resistance, the non- 
spherical (pear shape) shape of the earth, and the influence of other heavenly 
bodies are negligible. The law as applied to statcllites is as follows: The orbit of a 
satellite is an ellipse with the center of the earth at one focus. The path of a 
ballistic missile, not including the powered and reentry portions, is also an 
ellipse, but one that happens to intersect the surface of the earth. 

Kepler’s second law, or law of areas, states: Every planet revolves so that the line 
joining it to the center of the sun sweeps over equal areas in equal times. 

To fit earth orbital systems, the law should be restated thus: Every satellite 
orbits so that the line joining it with the center of the earth sweeps over equal 
areas in equal time intervals. 

When the orbit is circular, the application of Kepler’s second law is clear, as 
shown in Figure 5. In making one complete revolution in a circular orbit, a 
satellite at a constant distance from the center of the earth (radius r) would, 
for example, sweep out eight equal areas in the total time period (P = 1). Each 
of these eight areas is equal and symmetrical. According to Kepler’s second law, 
the time required to sweep out each of the eight areas is the same. When a 
satellite is traversing a circular orbit, therefore, its speed is constant. 

When the orbit is elliptical rather than circular, the application of Kepler’s 
second law is not quite so easy to sec; although the areas are equal, they arc 
not symmetrical (Fig. 6). Note, for example, that the arc of Sector I is much 
longer than the arc of Sector V. Therefore, since the radius vector sweeps equal 

•The terms “trajectory" and “orbit" are sometimes used interchangeably. Use of the term “trajectory’ 
came to astronautics from ballistics* the science of the motion of projectiles shot from artillery or firearms, 
or of bombs dropped from aircraft. The term ’’orbit" is used in referring to natural bodies, spacecraft, and 
manmade satellites. It is the path made by a body in its revolution about another body, as by a planet about 
the sun or by an artificial satellite about the earth. 
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Figure 5. Law of areas as applied to a circular orbit. 

areas in equal fractions of the total time period, the satellite must travel much 
faster around Sector T (near perigee) than around Sector V (near apogee). The 
perigee (a word derived from the Greek prefix peri-, meaning “near,” and the 
Greek root ge, meaning “pertaining to the earth”) is the point of the orbit nearest 

the earth. The apogee is that point in the orbit at the greatest distance from the / 

earth (the Greek prefix apo- means “from” or “away from”). 




Kepler’s third law, also known as the harmonic law, states: The squares of 
the sidereal periods * of any two planets are to each other as the cubes of their 
mean distances from the center of the sun. 

To fit an earth orbital system, Kepler’s third law should be restated thus: 
The squares of the periods of the orbits of two satellites are proportional to each 
other as the cubes of their mean distances from the center of the earth. The 



* The period of a planet about the sun. 
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moan distance is the length of the semimajor axis (a) of the ellipse, which is an 
average of the distances to perigee and apogee. Of course, in a circular orbit, 
the mean distance is the radius, r. 



Newton’s Laws 



While Kepler was working out his three laws of planetary motion, Galileo, an 
Italian physicist and astronomer, was studying the effects of gravity on falling 
bodies. Newton drew upon the work of both Kepler and Galileo to formulate 
his laws of motion. 

Newton’s first law states: Every body continues in a state of rest or of uniform 
motion in a straight line, unless it is compelled to change that state by a force 
imposed upon it. In other words, a body at rest tends to remain at rest, and a 
body in motion tends to remain in motion unless it is acted upon by an outside 
force. This law is sometimes referred to as the law of inertia. 

The second law of motion as stated by Newton says: When a force is applied 
to a body, the time rate of change of momentum is proportional to, and in the 
direction of, the applied force. If the mass remains constant, this law can be 
written as F = Ma. 

Newton’s third law of motion is the law of action and reaction: For every 
action there is a reaction that is equal in magnitude but opposite in direction to 

the action. If body A exerts a force on body B, then body B exerts an equal 

force in the opposite direction on body A. 

Force as Measured in the English System 

Newton’s three laws of motion are stated in terms of four quantities: force, 
mass, length, and time. Three of these, length, time and either force or mass, 
may be completely independent, and the fourth is defined in terms of the other 

three by Newton’s Second Law. Since the units and re. a* ; e values of these 

quantities were not known, Newton stated his second ’aw as a proportionality. 
Assuming that mass docs not change with time, this proportionality is stated as 
F x ma. )f proper units are selected, this statement may be written as an 
equation: 



F = ma 

The following are used in the metric system of measurement: 

F (dynes) = m (grams) times a (centimeters per second per second) 
F (Newtons) = m (kilograms) times a (meters per second per second) 



The most common force experienced is that of weight, the measure of the 
body’s gravitational attraction to the earth or other spatial body. Since this 
attraction is toward the center of the earth, weight, like any force, is a vector 
quantity. When the only force concerned is weight, the resulting acceleration is 
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normally called "g,” the acceleration due to gravity. For this special case, 
Newton's Second Law can then be written: 



W = Mg 



This equation can he used as a definition of mass. The value of g near the 
surface of the earth is approximately 32.2 feet per second per second: “g” is a 
vector quantity since it is directed always toward the center of the earth. If the 
weight, W, is expressed in pounds, rearranging gives: 

vr _ W (pounds) 

— g (feet/sec 2 ) 

The unit of mass in this equation is called a “slug.” Note that mass is a scalar 
quantity* and is an inherent property of the amount of matter in a body. Mass 
is independent of the gravitational field, whereas weight is dependent upon the 
field, the position in the field, and the mass of the body being weighed. 

Finally, Newton’s Second Law may now be written: 



F (pounds) = M (slugs) times a (ft/sec 2 ) 

The following example shows the use of this system of units and the magnitude 
of the “slug”: 



A package on earth weighs 161 pounds. 

Find: (a) its mass in slugs. 

(b) the force necessary to just lift it vertically from a surface. 

(c) the force necessary to accelerate it 10 ft/sec 2 on a smooth, level 
surface. 

(d) its weight if it were on the moon; assume the value of “g” there is 
% of that value here on the earth. 

(e) its mass on the moon. 

Solution: 

ra) M (slugs) = f- = = 5 siugs 

(D) F = W = M g = (5 slugs) (32.2 ft/sec 2 ) = 161 pounds 

The force must be applied upwards, in the direction opposite to weight. 

(c) F = Ma= (5 slugs) (10 ft/sec 2 ) = 50 pounds 

(d) W m „on = M g moon = (5 slugs) ( —g— ft/sec 2 ) = 26.83 pounds 

(e) M,noon “ “ 32T2 - ftTsec 2 ”” = 5 slu S s 



* A scalar quantity has magnitude only, in contrast to a vector quantity which has magnitude and direction. 
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This laat solution is, of course, to reemphasize that mass is independent of 
position. It will be shown later that the local value of g varies with altitude 
above the earth. It is significant to note that the weight will vary such that 
W 

the ratio remains constant. 

g 



Energy and Work 

Work, w, is defined as the product of the component of force in the direction of 
motion and the distance moved. Thus, if a force, F, is applied and a body moves a 
distance, s, in the direction the force is applied, w = Fs (Fig. 7). The units of work 
are foot-pounds. Work is a scalar as distinguished from a vector quantity. 



w - Fs 



r 

* 

i 
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Figure 7. Work performed as a force (F) is moved over the distance s. 



To do work against gravity, a force must be applied to overcome the weight, 
which is the force caused by gravitational acceleration, g. 

Therefore, F = Mg. If the body is lifted a height h (Fig. 8) and friction is negligi- 
ble, w = Mgh. For problems in which h is mich less than the radius from the cen- 
ter of the earth (h< <r), g may be considered a constant. 



M 



M 

7/,V/ ////*/// //////■ 



w=Fs 

F - Mg, s = h 
w = Mgh 



Figure 8. Work performed in lifting. 



If an object is pushed up a frictionless inclined plane, the work done is still Mgh 
(Fig. 9). 




Figure 9. Work performed on a frictionless inclined plane. 

For orbital mechanics problems, g varies and must be replaced by the value 
VgigT where the subscripts indicate the beginning and final values of g. In such 
cases 



w = MVgigs h 



Another type of work is that work done against inertia. If, in moving from one 
point to another, the velocity of a body is changed, work is done. This work against 
inertia is computed in the following steps: 



w 

but, F 
and 2as 

s 

so, w 



Fs 

Ma 

Vf 2 ~ V 0 2 

vr - v„ 2 
2a 

Fs = M a ill ~ v q 2 ) = M (yl - v,,-) = Myr _ Mv„- 

2a 2 2 2 



% X <) 

The quantity is defined as kinetic energy (KE). Therefore, work done 

against inertia (if the altitude and the mass remain the same) is equal to the change 
in kinetic energy. Energy is defined as the ability to do work, and it is obvious that a 
moving body has the ability to cU. work (for example, a moving hammer’s ability to 
drive a nail). A body is also able to do work because of its position or altitude; this 
is known as potential energy (PE). Units used to measure energy are similar to 
those used to measure work in that both are scalar rather than vector quantities. 

The sum of the kinetic and the potential energy of a body is its total mechanical 
energy. 
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Newton’s Law of Universal Gravitation 



Newton published his Principia in 1687 and included in it the law of universal 
gravitation, which he had been considering for about twenty years. This law 
was based on observations made bv Newton. Later work showed that it was 
onlv an approximation, but an extremely good approximation. The law states: 
Every particle in the universe attracts every other particle with a force that is 
proportional to the product of the masses and inversely proportional to the square 
of the distance between the particles. A constant of proportionality, G, termed the 
Universal Gravitational Constant, was introduced, and the law was written in this 
manner: 

r? _ Gmi m 2 

F 7 - — 



The value of G, the Universal Gravitational Constant, was first determined by Ca- 
vandish in a classical experiment using a torsion balance. The value of G is quite 
small (G = 6.6695 X 10 -s cgs units). In most problems the mass of one of the 
bodies is quite large. It is convenient, therefore, to combine G and the large mass, 
m,, into a new constant, p (mu), which is defined as the gravitational parameter. 
This parameter has different values depending upon the value of the large mass, mi. 
If m, refers to the earth, the gravitational parameter, p. will apply to all earth satel- 
lite problems. However, if the problem concerns satellites of the sun or other large 
bodies, p will have a different value based on the mass of that body. 

If we now simplify the law of gravitational attraction by combining G and mi and 
by adjusting the results for the English engineering unit system, we obtain the fol- 
lowing: 



G 




ft 3 

sec* 

m (Where F is lb force and m is slugs) 



If this expression is equated to the expression of Newton’s Second Law of 
Motion, as it applies in a gravitational field, we see that: 



F = mg 




m 



and after dividing by the unit mass, m, we obtain: 



g = 






r 2 



Thus, the value of g varies inversely as the square of the distance from the 
center o', the attracting body. 

For problems involving earth satellites, the following two constants are nec- 
essary for a proper solution: 



G nienrth 
r r (radius of earth) 



/W t h = 14.08 X 10** — 
20.9 X 10° ft 
2-15 



38 

ERIC 



The formulas must be used with proper concern for the units involved, and 
the value given for p. applies only to bodies attracted to the earth. 

Before applying Newton’s law of universal gravitation to the solution of prob- 
lems, it would be well to consider the possible paths that a body in unpowncd 
flight must follow through space. 

CONIC SECTIONS 

The conic sections were studied by the Greek mathematicians, and a body of 
knowledge has accumulated concerning them. They have assumed new signifi- 
cance in the field of astronautics because any jree-jlight trajectory can be 



Hyperbola 





Figure 10. Conic sections. 

represented by a conic section. The study of conic sections, or conics, is part of 
analytic geometry, a branch of mathematics that brings together concepts from 
algebra, geometry, and trigonometry. 

A conic section is a curve formed when a plane cuts through a right circular 
cone at any point except at the vertex, or center. If the plane cuts both sides 
of one nappe of the cone, the section is an ellipse (Fig. 10). The circle is a special 
case of the ellipse occurring when the plane cuts the cone perpendicularly to the 
axis. If the plane cuts the cone in such a way that it is parallel to one of the 
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sides of the cone, the section is called a parabola. If the phi'e cuts both 
nappes of the cone, the section is a hyperbola which has two branches. 

In one mathematical sense, all conic sections can be defined in terms of 
eccentricity (e). The numerical value of € is an indication of the relative shape 
of the conic (rotund or slender) and also an indication of the identity of the 
conic. 

If the eccentricity is zero, the conic is a circle; if the eccentricity is greater 
than zero but less than one, the conic is an ellipse; if the eccentricity is equal 
to one, the conic is a parabola; and if the eccentricity is greater than one, the 
conic is a hyperbola. 

Conic Sections and the Coordinate Systems 

In locating orbits or trajectories in space, it is possible to make use of either rec- 
tangular (sometimes called Cartesian) or polar coordinates. In dealing with artificial 
satellites, it is often more convenient to use polar rather than rectangular coordinates 
because the center of the earth can be used both as the origin of the coordinates and 
as one of the foci of the ellipse. 




Figure 11. Rectangular and polar coordinates superimposed on the conic sections. 



If rectangular and polar coordinates are superimposed upon a set of conics as 
shown in Figure 1 1 . equations of the curves can be derived. 
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The formula for the eccentricity c f a conic is e = Jj-. This ratio is constant for a 

specific curve. 

In Cartesian coordinates: 



r _ r _ V* 2 + y 2 
d k - x 

V x 2 + y- = e(k — x) 

Squaring both sides gives: 



x 2 + y 2 = e 2 (k — x) 2 



To convert the rectangular coordinates to polar coordinates, substitute as follows: 



x = r cos v (v is the lower case Greek nu) 
y = r sin v 

_ r r 

d k — r cos v 



ke — re cos v = r 
r + re cos v = ke 



r 



ke 

1 + e cos v 



This result is the general equation for all conics. 



Ellipse 

The ellipse is the curve traced by a point (P) moving in a plane such that the sum 
of its distances from two fixed points (foci) is constant. In the ellipse in Figure 12, 
the following are shown: the foci (F and F'); c, distance from origin to either focus; 
a, distance from origin to either vertex (semimajor axis) ; 2a, major axis; b, distance 
from origin to intercept on y-axis (suniminor axis); 2b, minor axis; and r 4- r', dis- 
tances from any point (P) on the ellipse to the respective foci (F and F r ). 

A number of relationships which are very useful in astronautics are derived 
from the geometry of the ellipse: 

r + r' = 2a (at any point on the ellipse) 
a 2 = b 2 + c 2 or a = \/b 2 + c 2 
b = \/a 2 — c 2 
c = \/a 2 - b 2 

c 

The eccentricity of the ellipse (e) = — .A chord through either focus perpendicu- 

2b 2 

lar to the major axis is called the latus rectum and its length = . 

3. 
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These relationships can be used to determine the parameters of an elliptical 
orbit of a satellite when only the radius of perigee and the radius of apogee 
are known. These parameters are important because, as is shown later, ey 



v 




Figure 1 2. Ellipse with center at origin of rectangular coordinate system. 

• 

are related to the total mechanical energy and total angular momentum of the 
satellite. Thereby they offer a means of determining these values through .the 
simple arithmetic of an ellipse rather than the vector calculus of celesttal 

mechanics. 

Sample problem: A satellite in a transfer orbit has a pengee at 300 NM 
above the surface of the earth and an apogee at 19,360 NM. Find a, b, c, and e 
for the ellipse traced out by this satellite. 

Solution: . 

Since the center of the earth is one focus of the ellipse, first convert the 
apogee and perigee to radii by adding the radius of the earth (3440 NM) . 



radius of perigee r p — 
radius of apogee r„ = 



altitude of perigee + radius of earth 
300 + 3440 = 3740 NM. 
altitude of apogee + radius of earth 
= 19,360 + 3440 = 22,800 NM. 
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Figure 13. Orbit of an artificial satellite showing radius of perigee and radius of 

apogee (not to scale). 



With this information, an exaggerated sketch of the ellipse can be made (Fig. 13). 
Compare this with Figure 12 to obtain: 



r p + r„ = major axis — 2a 
then 2a = 3740 + 22,800 = 26,540 NM 

or a = 26|40 _ 13>2 70 NM 

Also from comparing Figures 12 and 13: 

c = a - ip 

= 13,270 - 3740 = 9530NM 
Since a and c arc known, find b from the relationship given: 

b = V a ‘ ““ c 2 

or b = VCk 327 X 10 4 )- - (.95 3 X 10 4 )" 

b = V( 1-761 X 10 8 ) - (.908 X W) = V 853 X 10 s 
b = .923 X 10 4 = 9230 NM 



According to the formula given for eccentricity: 



c 



9530 

13270 



.718 



The ellipse is a conic section with eccentricity less than 1 (e < 1). 



Circle. The circle is a special case of an ellipse in which the foci have 
merged at the center; thus e = 0. The ellipse relationships can be used for a 
circle. 
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ENERGY AND MOMENTUM 



Once the basic geometry of a trajectory or orbit is understood, the next 
subject for investigation is the physics of energy and momentum. From concepts 
of linear and angular motion, concepts of linear and angular momenta log.caHy 
follow Once the formulas tor computing the specific angu ar 
the specific mechanical energy of a body in orbit arc delineated, then it a possible 
to solve for unknown quantities, such as the altitude of the y 
surface of the earth or the velocity at any point on the orbit. Any bo y P 
following a free-night path— whether it is a missile, a satellite, or a natural 
body— is governed by the laws of the conservation of specific mechanical energy 
and 5 ^ specific angular momentum. Once the value of either of these items is 
knowif at any point along a free-flight trajectory or orbit, then its value is known 
at all other points, since the value does not change unless the body is acted upon 

by some outside force. 



Mechanical Energy 



The law of Conservation of Energy states that energy can neither be created 
nor destroyed but only converted from one form to another. This law can be 
applied to orbital mechanics and restated in this way: The total mechanical 

energy of an object in free motion is constant, provided that no externa wor 
is done on or by the system. During reentry, work is done by the system and 
some of the mechanical energy is converted to heat. Similarly, during launch, 
work is done on the system as the propulsion units give up chemical energy. In 
this chapter, only the free-flight portion of the trajectory is considered, and it is 
assumed that there is no thrust and no drag. 

In order to establish a common understanding about changes in the amount 
of energy, it is necessary to agree upon a zero reference point for energy. 
Potential energy, or energy due to position, can be, and often is, measured from 
sea level. In working with earth-orbiting systems, however, the convention is 
consider a body as having zero potential energy if it is at an 
from the earth and as having zero kinetic energy if it is absolutely at rest v 
respect to the center of the earth. Under these circumstances, the total mechanical 
energy (PE + KE) is also equal to zero. If the total mechanical energy is 
positive — that is, larger than zero-the body has enough energy to escape from 
the earth. If the total mechanical energy is negative— that is, less than 
the body does not have enough energy to escape from the earth, and it mus 
be either in orbit or on a ballistic trajectory. 

The formula for PE, with the reference system as stated above, is PE 
-him . instead of using PE, a specific PE (PE per unit mass) can be used if both 

sides are divided by m; for example 



PE 

m 
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0 . 



If a body is at infinity, it has a specific PE equal to — ~£ ~ = „JL_ = 

r oo 

A similar case can be presented for kinetic energy. A body with some velocity 
relative to the center of the earth has kinetic energy defined by: 



KE = 



mv- 



Again, the specific kinetic energy (kinetic energy per unit mass) can be de- 
fined as: 



Specific KE = 



KE 

m 



v- 



In general, a body in free motion in space has a particular amount of me- 
chanical energy, and this amount is constant because of the conservation of 
mechanical energy. 

Total Mechanical Energy = KE + PE 

A more useful expression is obtained if we define Specific Mechanical Energy, 
E, or the Total Mechanical Energy per unit mass. Thus, we can write: 



E 

E 

E 



Total Mechanical Energy 

m 

KE PE 



m 

•> 

V“ 



m 



Specific Mechanical Energy, E, is also conserved in unpowered flight in space, 
ft- 

The units of E are — Since the mass term does not appear directly in the 

sec 

equation, E represents the specific mechanical energy of a body in general. 

If the solution to the Specific Mechanical Energy equation yields a negative 
value for E, the body is on an elliptical or circular path (nonescape path). If 
E is exactly equal to zero, the path is parabolic; this is the minimum energy 
escape path. If E is positive, the path is hyperbolic, and the body will also escape 
from the earth's gravitational field. 

Although the value of E, once determined, remains constant in free flight, there 
is a continuous change in the values of specific kinetic energy and specific po- 
tential energy. High velocities nearer the surface of the earth, representing high 
specific kinetic energies, are exchanged for greater specific potential energies 
as distance from the center of the earth increases. In general, velocity is traded 
for altitude, kinetic energy is traded for potential energy. The sum remains 
constant. 

Linear and Angular Momentum 

When a body is in motion, it has momentum. Momentum is the property a 
body possesses because of its mass and its velocity. In linear motion, momentum 

is expressed as mv and has the units, foot-slug, 

sec 
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When a rigid body, such as a flywheel, rotates about a center, it has angular 
momentum. Once the flywheel is in motion, its angular momentum would remain 
constant if it were not acted upon by forces such as friction and air resistance. 
Similarly, a gyroscope would rotate indefinitely in the absence of friction and air 
resistance. Thus, ignoring such losses, angular momentum will remain constant. 
In space, it can be assumed that such forces are negligible and that angular 
momentum is conserved. This is another tool to use in analyzing orbital systems. 

Angular momentum is the product of moment of inertia, I, and the angular 
velocity, co. Moment of inertia of a body of mass, m, rotating about a center 
at a distance, r, can be expressed as m r J . The angular momentum is then equal 

to m r- co. . 

For convenience in calculations, the term Specific Angular Momentum, H, is 

defined as the angular momentum per unit mass. Remembering that the magm- 



Angular Momentum - m r 2 o) 



horizontal 
f.o _ 



Ang. Momentum 



,*.H = v u r (Circular Motion) 
h 

Figure 14. Specific angular momentum of a circular orbit. 

tude of the instantaneous velocity vcoor of a body rotating in constant circular 
motion about a center with radius r is equal to cor and that the vector is perpen- 
dicular to the radius, the expression for specific angular momentum of a circular 
orbit can be simplified as shown in Figure 14. 

The general application of specific angular momentum to all orbits requires 
that the°component of velocity perpendicular to the radius vector be used. This 
velocity component is defined as 

V h = V COS (£> 

where (f> is the angle the velocity vector makes with the local horizontal, a line 
perpendicular to the radius. In an elliptical orbit, the geometry is as shown m 
Figure 15. The body in orbit has a total velocity v which is always tangent to the 
flight path. 

The formula, H = v r cos </>, defines the specific angular momentum for all or- 
bital cases. The angle <f> w ihz flight path angle and is the angle between the local 
horizontal and the total velocity vector. It should be noted that the angle <p is 
equal o zero for circular orbits. Further, in elliptical orbits, (f> is zero at the 
points of apogee and perigee. 

The two important formulas that have been presented in this section are those 
for E and H. These formulas permit a trajectory or an orbit to be completely 
defined from certain basic dat 
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V- 



ft- 

sec 2 



E 






r ’ 



when units of E are 



H = 



v r cos <f>, when units of H aie 



ft 2 

sec 



If v, r, and </> are known for a given trajectory (or orbit) at a given position, 
then E and H can be determined. In the absence of outside forces, E and H 
are constants; therefore v, r, and <f> can be determined at any other position on 
the trajectory or orbit. Equations for the specific angular momentum and the 
specific mechanical energy can be used in practical application to the two-body 
problem and to the free-flight portion of the ballistic missile trajectory. 




Figure 15 . Specific angular momentum. 



THE TWO-BODY PROBLEM 

It is implicit in Newton’s law of universal gravitation that every mass unit 
in the universe attracts, and is attracted by, every other mass unit in the universe. 
Clearly, small masses at large distances are infinitesimally attracted to each other. 
It is neither feasible nor necessary to consider mutual attractions of a large 
number of bodies in many astronautics problems. The most frequent problems of 
astronautics involve only two interacting bodies: a missile payload, or satellite, 
and the earth. In these instances, the sun and moon effects are negligible except 
in the case of a space probe, which will be noticeably affected by the moon, if 
it passes close to the moon, and which will be controlled by the sun, if it escapes 
from the earth’s gravitational field. 

Military officers concerned with operational matters are primarily interested 
in launching a missile from one point on the earth’s surface to strike another point 
on the earth’s surface and in launching earth satellites. In these problems, the 
path followed by the payload is adequately described by considering only two 
bodies, the earth and the payload. The problem of two bodies is termed the 
two-body problem; its solution dates back to Newton. 

It is indeed fortunate that the solution of the two-body trajectory is simple 
and straightforward. A general solution to a trajectory involving more than two 
bodies does not exist. Special solutions for these more complex trajectories 
usually require machine calculation. 
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Figure 16. Trajectory relationships. 

The two-body problem is described graphically in Figure 16. A small body, 
nv, has a velocity, v, at a distance r from the origin chosen as the center of 
mass of a very. massive body, nil. The problem is to establish the path followed tty 
body, nr., or to define its trajectory. This is a typical problem in mechanics 
given the present conditions of a body, what will these conditions be at any 
time, t. later? First, we shall find r as a function of v, where v is the polar 
angle measured from a reference axis to the radius vector. 

At the outset, it should be apparent that the entire trajectory will take place 
in the plane defined by the velocity vector aRd the point origin. There are no 
forces causing the body, nr, to move out of this plane; otherwise, the conditions 
arc not those of a two-body free-flight problem. 

In the earlier outline of the laws of conservation of energy and momentum, 
the following conditions were established: 

- V- = E = a constant (D 

2 r 

H = vr cos (J> = a constant (2) 

Equations (1) and (2) can be combined and, with the aid of the calculus, 
the following equation can be derived: 



r 



H 2 /M 



1 + 




2EFF 

M" 



cos V 



(3) 



Equation (3) is the equation of a two-body trajectory in polar coordinates. 

Earlier, the following equation was given as the equation of any conic section in 
polar coordinates, the origin located at a focus: 

_ ke 

r — 1 + e cos v (4) 

Equations (3) and (4) are of the same form; hence, equation (3) is also the 
equation of any conic section (origin at a focus) in terms of the physical constants, 
E and H, and the two-body trajectories are then conic sections. This conclusion sub- 
stantiates Kepler’s first law. In fact, Kepler’s first law is a special case because an 
ellipse is just one form of conic section. 
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Since equations (3) and (4) are of the same form, it is possible to equate like 
terms, which will lead to relationships between the physical constants, E, H, and /x, 
and the geometrical constants, e, a, b, and c. Thus: 



and 




€ 




2EH- 



(5) 



( 6 ) 



Physical Interpretation of the Two-Body Trajectory Equation 

Analysis of the two-body trajectory equation will give an understanding of the 
physical reaction of a vehicle (small body) under the influence of a planet 
(large body). 

If E < 0, the trajectory is an ellipse What is the condition that E be less than 
zero? It is simply that the kinetic energy of the small mass, m^>, because of its 
relatively low velocity, is less than the magnitude of its potential energy. There- 
fore, the body cannot possibly go all the way to infinity; that is, it cannot go 
to a point where it is no longer attracted by the larger body — where the potential 
energy is zero. The smaller mass cannot escape. It must remain “captured” by 
the force field of the larger body. Therefore, it will be turned back toward the 
larger body, or, more in keeping with the idea of potential energy, it will always 
be ‘falling back” toward the more massive body. When this particular balance 
of energy exists, the trajectory is elliptical with one focus coincident with the 
center of mass of the larger body. In the actual physical case, the larger body 
will have a finite size; that is, it will not be a point mass, and this ellipse may 
intersect the surface of the larger body as it does in the case of a ballistic missile. 
If the velocity is sufficiently high, and its direction proper, the ellipse may com- 
pletely encircle the central body, the condition of a satellite. 

If E = 0, the kinetic energy exactly equals the magnitude of the potential 
energy, and the small mass. nv>. has just enough energy to travel to infinity, away 
from the influence of the central body, and come to rest there. The small body 
will follow a parabolic path to infinity. The velocity which is associated with this 
very special energy level is also very special and is commonly called the “escape 
velocity.” 

Escape velocity can be calculated by setting E = 0 in the mechanical energy 
equation (1) as follows: 



V^csp 

2 




= E = 0 



esc 




V 



esc 
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Thus, it can he seen that escape velocity decreases with distance from the 



center of the earth. At the earth’s surface, 




(2) (14.08) (10 ir> - 
(20.9) ( 10* ft) 




v esr = 36,700 ft/sec. 

If the velocity of the small mass exceeds escape velocity, which will be the 



would still have velocity at infinity. In a physical sense, it would have velocity 
relative to the earth at any large distance from the earth. 

Considering the sounding rocket, only the straight-line, degenerate conic is a 
possible trajectory. But. again, the value of E will determine whether escape is 
possible; that is, if E < O, the straight-line trajectory cannot extend to infinity. 
If E = O or E > O. the straight line will extend to infinity. 

Example Problem'. The first U.S. * l moon shot, the Pioneer I, attained a height 
of approximately 61.410 NM above the earth's surface. Assuming that the 
Pioneer had been , mounding rocket (a rocket fired vertically), and assuming 
a spherical, nonrotating earth without atmosphere, calculate the following: 

a. E (total specific energy) 

b. Impact velocity (earth’s surface) 

Solution: Given 

(a) At apogee (greatest distance from earth) : 

Altitude (above earth's surface) = 61,410 NM 

Earth radius = 3440 NM 

Velocity = O (Only for a sounding rocket) 

r = altitude 4- earth’s radius 

r = (61,410 + 3440) NM = 64,850 NM 




the mass will have some residual velocity. In a mathematical sense, the body 



v- 



(14.08) (10 15 ) ~ 



E 



2 



0 



r 



(64,850) NM (6080) 



E = - 3.57 X 10 7 ft-/sec- 



Answer 



(b) Since the specific energy is constant, 
At the earth’s surface: 



r = 3440 NM 
E = - 3.57 X 10 7 ftVsec- 
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E 



it 

r 



v~ 

"2 



3.57 x 10 T 



ft- 

SCC' 



V- 



(14.08) (lO 15 -^--) 



sec- 



(3440) NM (6080) 



ft 



v- = 2 f ( — 3.57 X 1 0 7 ) + ( 67.4 X 10“)] = 2(63.8 X 10 T )ft a /scc 2 

v limi.m-n = 35.700 ft/sec Answer 



This is also the approximate burnout velocity of the vehicle. As the surface 
escape velocity is 36.700 ft/sec, it is clear that Pioneer I did not attain escape 
velocity, and so it returned to earth. 



Elliptical Trajectory Parameters 

While parabolic and hyperbolic trajectories, especially the latter, are of interest in 
problems of interplanetary travel, elliptical trajectories comprise the ballistic missile 
and satellite cases, which are of current military interest. It is important, then, to re- 
late the dimensions of an ellipse (a, b. and c) to the physical constants (E, H, and 
/x) as was previously done for e. 

The relationship, r (1 + r,, = 2a, was presented earlier. If this equation is applied 
to point P in Figure 17, 




Figure 17. Ellipse. 



r = — — — and r p = a — c. Then, 
1 + e 



ke 

1 + e 



a — c 



(S) 
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But e = ; therefore, substituting c = ea into (8), 



— = a — ea = a(l — e) 
1 + e 

ke = a(l — €-) 

But from equations (5) and (6), 



( 9 ) 



ke = W/ix and € 



-V 



1 + 



2EH 2 

•> 

/x- 



Substituting these relationships into equation (9), 

HI = a (l - 1 - m~) 

i V /*- / 



H 2 2EH 2 a 



_ 1 = 2Ea 






a 



_ _ A 



2E 



(EXTREMELY USEFUL) 



Also a = 



From the following equation: 






’(4 - +) 



( 10 ) 

( 11 ) 



Jbi 

a 2 

_bi 

a 



= 1 — e 2 



= a(l — e 2 ) 



But from equation (9), 



a(l - e 2 ) = ke = 



H 2 



Therefore, 



b 2 . H 2 



a 



( 12 ) 



Equations (10) and (12) are extremely important relationships; an understand- 
ing of them is essential to material that follows on ballistic missiles and satellites. If 
injection conditions of speed and radius are fixed, it is clear that a, the semimajor 
axis of the elliptical trajectory, becomes fixed, regardless of the va^ae of the flight 
path angle at burnout. Equation (10) points out there k a direct relationship be- 
tween the size of an orbit and the energy level of the orbiting object. . juation (12) 
points out that for a given energy level there is direct relauonshm between the 
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length of the semi latus rectum of an elliptical trajectory (a shape parameter) and 
the specific angular momentum of the orbiting object. This implies that the size and 
shape of an elliptical trajectory are determined by the E and H. 



Two-Body Trajectory Definitions and Geometry 

The general equation of two-body trajectories has now been introduced. Before 
proceeding to problem applications it would be well to consider in detail some 
commonly used terms and symbols. 

First, refer to Figure 18. In general, the point P is called the periapsis and 
P' the apoapsis. If the earth is at point 0, the ellipse would then represent the 
trajectory of an earth satellite; P is then termed perigee and P' apogee. If 
the sun is at point O, the ellipse would represent a planetary orbit; P is then 
called perihelion and P' aphelion. 

In order to explain the use of the angle v, the geometry of satellites will be 
discussed briefly. Figure 18 depicts a planetary orbit (not to scale). 

In astronomy and celestial mechanics it is standard practice to measure a body’s 
position from perihelion point P. There are several reason.^ for using pctihelion, 
including the fact that perihelion of any body in the y-Air system excop 1 M"’c:u r v 
and Venus is closer to the earth’s orbit than is the i' <.‘ '■ nj IvT/ a in 
fm - a highly ccccnti"' orbit such as a comet’s, the body would not b- vUiMe 
at aphelion, in order to conform to accepted practice, tlv t, the angle nu (v), 
measured fr? n per! p : ..s, has been introduced. This angle, which is called the 
true anomaly, is of considerable importance in time-of-flight calculations. 




The true anomaly, however, does not lend itself well to ballistic missile problems, 
as can be seen from the simplified ballistic missile geometry in Figure 19. 
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Figure K -JUstSc missile trajectory. 



In a ballistic missile trajectory, perigee is entirely fictitious. The missile obviously 
never traverses the dashed portion of the trajectory. The solid portion of the trajec- 
tory is all that is of real interest, and this portion is in the second and third quad- 
ram.'. of the angle v. It is convenient then to define an angle 9, measured counter- 
clockwise from apogee (apoapsis, in general), such that 

v = 0 + it. (13) 

6 will normally have values in the first and fourth quadrants. From ( 13) , it is clear 
that derivatives of v and 9 will be interchangeable. The equation of a conic section 
in terms of 9 can be found by substituting (13) into (4), 

_ ke = kf 

r l+e cos v 1 + e cos (6 4- tt) 

r = . 

1—e cos 6 

With this understanding of the relationship between v and 0 , it will be convenient 
to use v when working with satellite and space trajectories and 6 when working with 
ballistic missiles (see App. C). 



EARTH SATELLITES 

During their free flight trajectory, satellites and ballistic missiles follow paths 
described by the two-body equation. For a satellite to achieve orbit, enough 
energy must be added to the vehicle so that the ellipse does not intersect the 
surface of the earth. However, not enough energy is added to allow the vehicle 
to escape. Therefore, the ellipse and the circle are the paths of primary interest. 

The orientation, shape, and size of orbits are important to the accomplishment 
of prescribed missions. Therefore, eccentricity (e), major axis (2a), minor axis 
(2b), and distance between the foci (2c) are of interest. It is necessary to know 
the relationships of these geometric values to the orbital parameters in order to 
make a.i analysis of an orbit. For example, it is helpful to remember that: 
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r,, (radius at perigee) — a — c 
r., (radius at apogee) = a *f c 
r r: + r n = 2a 

e = 

a 

a- = b- + c- 



Specific mechanical energy. E. and specific angular momentum, H, are of primary 
concern when elliptical and circular orbits are discussed. If there are no outside 
forces acting on a vehicle in an orbit, the specific mechanical energy and the specific 
angular momentum will have constant values, regardless of position in the orbit. 
This means that if E and H are known at one point in the orbit, they are then known 
at each and every other point in the orbit. At a given position if radius r, speed v, 
and flight path angle $ are known, E and H can be determined from: 



P _ jv^_ fX_ _ fX 

n ■” 2 r 2a 

H — vr cos <j> 

If the values of E and H are known for a particular orbit, and the speed and 
flight path angle at a certain point in the orbit are to be determined, the energy 
equation can be solved for v, and then the angular momentum equation can 
be solved for 0. 

The equations for the speed in circular and elliptical orbits are important. 
The equation for circular speed is: 




The equation for elliptical speed is: 



v = 




Another equation that is important in the analysis of orbits is the equation 
for orbital period. For a circular orbit the distance around is the circumference 
of the circle which is 2m. Therefore, the period, which is equal to the distance 
around divided by the speed, is this: 

P = 2?rr 

v * 

Now, substitute for v the speed in circular orbit: 
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V 



= V-f- 

2nT 



i M 
\ r 



Multiply numerator and denominator of the right hand side by V r: 



P = 



Vr 



yJi 



2wr“*' : 

Vm 



Using the principle of Kepler’s third law, replace r by the mean distance from 
the focus, which is equal to the semimajor axis a, and the equation becomes: 



P = 



27ra : i; - 

Vm 



Squaring both sides. P- = . Since is a constant, P- is proportional 

M ( sec- \ 

to a 3 , and for earth satellites P- = { 2.805 X lO" 1 " J (a) 3 . Or, P = 

5.30 X 10-*- s ?A (a ){ 
ft- / 

Problem: Initial data from Friendship 7 indicated that the booster burned out 
at a perigee altitude of 100 statute miles, speed of 25,700 ft/sec, and flight path 
angle of O'. Determine the speed and height at apogee, and the period. 

Given*: h p = 100SM = .5 X 1 0 R ft r e = 20.9 X 10° ft 

v b0 = 25.700 ft/sec 

</>bo = 0° 

Find: v u , h„. P 





Figure 20. Orbit of Friendship 7 (not to scale). 



* Even though <£ bo = 0°, if burnout altitude were not given as perigee altitude, you would have to de- 
termine if this were perigee or apogee. To do this, you would compute the circular speed for the given 
burnout altitude and compare thus with the actual speed. If the circular speed is greater than the actual, 
burnout was at apogee; if the circular speed was less than the actual speed, burnout was at perigee. 
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Solution: 



E = 



v- 

2 



P- 



(2.57 X 10»)- 14.08 X 10 15 

2 (20.9 + .5) X 10 6 



= (3.31 X 10 8 ) - (6.58 X 10 8 ) = -3.27 X 10 8 



but E = - 2a = - 



P _ 



-14.08 X 10 1S 
-3,27 X 10 8 

= 43.1 X 10® 



From Figure 20, r n + r p = 2a 

.'.r a = 2a - r p « (43.1 - 21.4) X 10 8 
= 21.7 X 10 6 

h„ = r a - r c = (li.7 - 20.9) X 10 8 = .8 X 10 8 fl 

= 151 sm Answer 



H p = H a 



v p r p - v n r„ 
V a = 



v„ = 



_ v p r P 



p Ap y 

(2.57 X 10 *) (21.4 X 10 6 ) 



21.7 X 10° 



25,400 ft/sec 



A 4> H S6C 

P- = _ 47r ~ a — = (2 805 X 10~ 15 ) (21.6 X 10°ft) 3 -|^ 



= 17,300 mph Answer 
sec 2 
ft 3 

= 28.1 x 10 6 sec 2 



P = 5.30 X 10 3 sec = 88.3 min Ans. 

It is interesting to compare the computed apogee and period results with the 
actual orbit (later data gave a higher accuracy for burnout conditions): 



Item Actual figures Computed figures 

v b0 25,728 ft/sec 25,700 ft/sec 

h b0 97.695 SM 10CSM 

h. 158.85 SM 151 SM 

p 88.483 min 88.3 min 



Note that using three significant figures results in h p = .5 X 10° ft, about 94 
SM. Such errors are common using slide rule accuracy, but this problem does 
illustrate the techniques used. 

From the foregoing problem, it is evident that the principles and relatively 
simple algebraic expressions presented thus far are extremely important. *hey 
enable one a) analyze a trajectory or orbit rather completely — with a slide rule 
for academic or generalized discussion purposes, or with a digital computer for 
system design and operation. The discussion has been, however, confined to the 
two dimensional orbital plane. Before discussing some of the more interesting 
facets of orbital mechanics, it is necessary to properly locate a payload in three 
dimensions. 
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LOCATING BODIES IN SPACE 



Tn one of the coordinate systems for space used by engineers and scientists, 
the origin is the center of the earth. This is a logical choice since the center of 
the earth is a focus for all earth orbits. 

With the center of the coordinate systems established, a reference frame is 
required on which angular measurements can be made with aspect to the center. 
The reference frame should be regular in shape, and it should be fixed in space. 
A sphere satisfies the requirement of a regular shape. The sphere of the earth 
would be a handy reference if it were fixed in space, but it rotates constantly. 




Figure 21. Celestial Sphere. 




Figure 22. The vernal equinox. 



1 herefore. the celestial sphere is used to satisfy the requirement for a reference 
frame. This is a nonrotating sphere of infinite radius whose center coincides with 
the center of the earth and whose surface contains the projection of the celestial 
bodies as they appear in the sky (Fig. 21). The celestial equator is a projection 
ol the earth's equator on the celestial sphere. The track of a satellite can be 
projected on the celestial sphere by extending the plane of the oribt to its inter- 
section with the celestial sphere. 

After the center of the system and the celestial equator have been defined, 
a reference is required as a starting point for position measurements. This point, 
determined at the instant winter changes into spring, is found by passing a line 
from the center of the earth through the center of the sun to the celestial equator, 
and is called the vernal equinox (Fig. 22). 

After the references for the coordinate system have been established, the orbit 
itself must be located. The first item of importance is right ascension (f!) of the 
ascending node, which is defined as the arc of the celestial equator measured 
eastward from the vernal equinox to the ascending node (Fig. 21). The ascend- 
ing node is the point where the projection of the satellite path crosses the celestial 
equator from south to north. In other words, right ascension of the ascending 
node is the angle measured eastward from the first point of Aries to the point 
where the satellite crosses the equator from south to north. 

The next item ot importance is the angle the path of the orbit makes with the 
equator. This is the angle of inclination (i), which is defined as the angle that 
the plane of the orbit makes with the plane of the equator, measured counter- 
clockwise from the equator at the ascending node. Equatorial orbits have i = 0°; 
posigrade orbits have i = 0° to 90 5 ; polar orbits have i = 90 c ; and retrograde 
orbits have i= 90° to 180°. 
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To describe the orbit further, the perigee is located. The angular measurement 
from the ascending node to the perigee, measured along the path of the orbit in 
the direction of motion, is called the argument of perigee ( co ). 

If, in addition to the coordinates of the orbit, a time of either perigee or right 
ascension of the ascending node is known, along with the eccentricity and the 
major axis of the orbit, the exact position and velocity of the satellite can be 
determined at any time. Six quantities (right ascension, inclination, argument of 
perigee, eccentricity, major axis, and epoch time at either perigee or ascending 
node) form a convenient grouping of the minimum information necessary to 
describe the orbital path as well as the position of a satellite at any time. They 
constitute one set of orbital elements, known as the Breakwell Set of Keplerian 
Elements. 



Orbital Plane 

Another interesting facet of earth satellites concerns the orbital plane. There 
is a relationship between the launch site and the possible orbital planes. This 
restriction arises from the fact that the center of the earth must be a focus of 
the orbit and, therefore, must lie in the orbital plane. 

The inclination of the orbital plane, i, to the equatorial plane is determined 
by the following formula: cos i (inclination) = cos (latitude) sin (azimuth) where 
the azimuth is the heading of the vehicle measured clockwise from true north. 

As an example, a satellite launched from Cape Kennedy and injected at 30° N on 
a heading due east (Azimuth 90°) will lie in an orbital plane which is inclined 30° 
to the equatorial plane. 



cos i “ (cos latitude) (sin azimuth) 
= cos 30° sin 90° 
cos i = cos 30° 
i = 30° 



It can be deduced from the above that the minimum orbital plane inclination for a 
direct (no dog leg or maneuvering) injection will be closely defined by the latitude 
of the launch site. All launch sites, thereby, will permit direct injections at inclina- 
tion angles from that minimum (the approximate latitude of the launch site) to 
polar orbits (plus retrograde supplements), provided there were no geographic re- 
strictions on launch azimuth, such as range safety limitations. For example, direct 
injections from Vandenberg AFB (35°N) would permit inclination angles from 
about 35° to 145°. 

Once the inclination of the orbital plane is defined, the ground track can be dis- 
cussed. 
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SATELLITE GROUND TRACKS 



The orbits of all satellites lie in planes which pass through the center of a 
theoretically spherical earth. Each plane intersects the surface of the earth in a 
great circle (Fig. 23). 




Figure 23. Satellite ground track geometry. 

A satellite’s ground track is formed ' y the intersection of the surface of the 
earth and a iine between th rt center of the earth and the satellite. As the space 
vehicle moves in its orbit, this intersection traces out a path on the ground below. 

There are five primary factors which affect the ground track of a satellite 
moving along a free flight trajectory. These are: 



1. Injection point 




7. Inclination angle 


(i) 


3. Period 


<P) 


4. Eccentricity 


(e) 


5. Argument of Perigee 


(«) 



Of the above, the injection point simply determines the point on the surface from 
which the ground track begins, following orbital injection of the satellite. Inclination 
angle has been discussed in the previous section and will be treated below in further 
d r *ail. Period, eccentricity, and argument of perigee each affect the ground track, 
bi..t it is often difficult to isolate the effect of any one of the three. Therefore, only 
general remarks regarding the three factors will be made, rather than an intricate 
mathematical treatment. 

If the study of satellite ground tracks is predicated upon a nonrotating earth, 
the track of a satellite in a circular orbit is easy to visualize. When the satellite’s 
orbit is in the equatorial plane, the ground track coincides with the equator 
(Fig. 24). 
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Figure 24. Equatorial track. 

If the plane of the orbit is inclined to the equatorial plane, the ground track 
moves north and south of the equator. It moves between the limits of latitude 
equal to the inclination of its oroital plane (Fig. 25). A satellite in either circular 
or elliptical orbit will trace out a path over the earth between these same limits 
of latitude, determined ny the inclination angle. However, the satellite in elliptical 
orbit will, with one exception, remain north or south of the equator for unequal 
periods of time, "Hiis exception occurs when the major or long axis of the orbit 
lies in the equatorial plane. 




Figure 25. North-South travel limits. 

The inclination of an orbit is determined by both the latitude of the vehicle and 
the direction of the vehicle’s velocity at the time of injection or entry into orbit. 
That is, the cosine of the inclination angle equals the cosine of the latitude times 
the sine of the azimuth (when the azimuth is measured from north). The mini- 
mum inclination which an orbital plane can be made to assume is the number of 
degrees of latitude ai which injection occurs. Tuis minimum inclination occurs 
when the direction of the vehicle’s velocity is due east or west at the time of in- 
jection. If the vehicle’s direction at injection into orbit is any direction other than 
east or west, the inclination of the orbital plane will be increased (Fig. 26). 
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Figure 26 . Injection-inclination geometry. 

On a flat .’nap of the earth, satellite ground tracks apnear to have different 
shapes than on a sphere. The ground track for a vehicle ? an inclined circular 
or elliptical orbit appears as a sinusoidal trace with N orth-Vcuth limits equal to 
the inclination of the orbital plane (Fig. 27). 




Figure 27. Ground track on flat, non-rotating earth map. 

When the earth’s rotation is considered, visualizing a satellite’s ground track 
becomes more complex. A point on the equator moves from west to east more 
rapidly than do points north and south of the equator. Their speeds are the 
speed of a point on the equator times the cosine of their latitude. Satellites in 
circular orbit travel at a constant speed. However, when the orbits are inclined 
to the equator, the component of satellite velocity which is effective in an easterly 
or westerly direction varies continously throughout the orbital trace (Fig. 28). 
As the satellite crosses the equator, its easterly or westerly component of velocity 
if its instantaneous total velocity times the cosine of its angle of inclination. 
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Figure 28. Effective East/West component of satellite velocity. 

When it is at the most northerly or southerly portion of its orbit, its easterly or 
westerly component is equal to its total instantaneous velocity. 

In elliptical orbits only the horizontal velocity component contributes to the 
satellite's ground track. Further complication results because the inertial or ab- 
solute speed of the satellite varies throughout the elliptical path (Fig. 29). 

Because the ground track is dependent upon the relative motion between the 
satellite and the earth, the visualization of ground tracks becomes quite compli- 
cated. Earth rotation causes each successive track of a satellite in a near earth 
orbit to cross the equator west of the preceding track (fig. 30). This westeoy 
regression is eaual to the period of the satellite times the rotational speed of the 
earth. The regression is more clearly seen if angular speed is considered. The 
earth’s angular speed of rotation is 15°/hour. The number of degrees of regres- 
sion (in terms of a shift in longitude) can be determined by multiplying the 
period of the satellite by 15°/hour, the angular speed of the earth. If the altitude 
of a satellite is increased, t’n reby increasing the time required to complete one 
rcvo.ution in the orbit, the distance between successive crossings of the equator 
increases. 




Figure 29. Variation of velocity magnitude in an elliptical orbit. 
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Figure 30. Ground track regression due to earth rotation. 



Again, using a flat map of the earth, the track of a satellite in circular orbit (with 
a period of less than 20 hrs) appears as a series of sinusoidal traces, each succes- 
sively displaced to the west (Fig. 31). 

The ground track of a satellite in elliptical orbit results in a series of irregular 
traces on a flat map which have one lobe larger than the other. The lobes are 
compressed by an amount which depends on orbital time, are altered in shape 
by the combined factors (inclination, eccentricity, period, and location of perigee), 
and successive traces are displaced to the west (Fig. 32). 

Some satellites follow orbits that have particularly interesting ground tracks. 
A satellite with a 24-hour period of revolution is one such case. If this satellite is 
in a circular orbit in the equatorial plane, it is often referred to as a synchronous 




Figure 31. Westward regression of sinusoidal tracks. 




satellite; its trace is a single point. If it orbits in the polar plane, it will complete 
half of its orbit while the earth is rotating halfway about its axis. The result is a 
trace which crosses a single point on the equator as the satellite crosses the 
equator heading north and south. The complete ground track forms a figure eight. 
If the p) \ne of the circular orbit is inclined at smaller angles to the equator, 
the figure eights are correspondingly smaller (Fig. 33). 




Figure 33 Figure eights for inclinations less than polar. 



The shape of the figure eighi may be altered by placing a satellite in an 
elliptical flight path. The eccentricity of the ellipse changes the relative size of 
the loops of the figure eight. If eccentricity and inclination are fixed, then changrng 
the location of perigee will vary the shape and orientation of the figure eight as 
shown in Fig. 34. The longitude of perigee is, of course, determined by the con- 
ditions and geographical location of injection into the 24-hour orbit. However, the 
latitude of perigee is fixed by the inclination of the orbital plane (i), and the argu- 
ment of perigee (cu) : 

Sin (latitude of perigee) = sin i sin 

Apogee is located on the same meridian as perigee and is at the same degree of 
latitude but in the opposite hemisphere. 

Circular ground tracks similar to the one shown in Fig. 34(g) have been 
proposed in certain navigation satellite concepts. In this case with apogee in the 
northern hemisphere, a multiple satellite system using this type of ground track 
might be used for supersonic aircraft navigation in the North Atlantic. 
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Eccentricity (r) =• .6 
Inclination (i) - 60° 

Argument of Perigee (w) = As Stated 



n're 34. Variation of elliptical, 24*hour track with movement of perigee. 
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L = 37.75° S 
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Eccentricity (f) * .6 
Incline rioi; (i) = 60° 

Argumen: of Perigee (<j) = As Stated 



Figure 34. Variation of elliptical, 24-hour track with movement of perigee, continued. 
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If a satellite is in an orbit with a period much greater than the earth’s twenty 
four hour period of rotation, the satellite appears as a point in space under which 
the earth rotates. If the orbit is in the equatorial plane, the trace is a line on the 
equator moving to the west. If the orbit is inclined to the equator, its *arth track 
will appear to be a continuous trace wound around the earth like a spiral between 
latitude limits equal to its angle of inclination (Fig. 35). An example of this 
phenomena is the ground track of the moon. 

It should now be clear that there is an almost limitless variety of satellite ground 
tracks. To obtain a particular track, it is only necessary that the proper orbit be 
selected. If changes are made in the inclination of an orbital plane to the 
equator, if its period is varied, if the eccentricity is controlled, or, if the location 
or perigee is specified, many different ground tracks can be achieved. 




/ 

Figure 35. Track for inclined orbit with period greater than one day. 



SPACE MANEUVERS 

One characteristic of satellites is that their orbits are basically stable in inertial 
space. This stability is often an advantage, but it can also pose problems. Space 
operations such as resupply, rendezvous, and interception may require that the 
orbits of space vehicles be changed. Such changes are usually changes in orbital 
altitude, orbital plane (inclination), or both. In this section some methods of 
maneuvering in space will be reviewed. 

Altitude Change 

When a satellite or space vehicle is to have its orbit changed in altitude, addi- 
tional energy is required. This is true whether the altitude is increased or de- 
creased. The classic example of changing the orbital altitude of a satellite is the 
HOHMANN TRANSFER. 

The Hohmann transfer is a two-impulse maneuver between two circular, 
coplanar orbits. For most practical problems, this' method uses the least amount 
of fuel and is known as a minimum energy transfer. The path of the transfer 
foliows an ellipse which is cotangential to the two circular orbits (Fig. 36). 
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In accomplishing a Hohniann transfer, two applications of thrust arc required. 
Each application ot thrust changes the speed of the vehicle and places it into a 
new orbit. Obviously both the direction and magnitude of the velocity change, Av, 



4v 




Transfer 

ellipse 



must be accurately controlled for a precise maneuver. If an increase in altitude is 
desired, the point of departure becomes the perigee of the transfer ellipse; the 
point of injection to the higher circular orbit becomes the apogee of the transfer 
ellipse. (For the transfer ellipse, 2a = ri + ro.) To lower altitude, the reverse 
is true. The point of departure will be the apogee of the transfer ellipse. 

In general, the process for determining the total increment of velocity, Av, 
required to complete a Hohmann transfer can be divided into seven steps. 



(1) Determine the velocity the vehicle has in the initial orbit. 



v - = \Hr 

(2) Determine the velocity required at the initial point in the transfer orbit. 



v. = JJt - -t. 

\ ri a 

(3) Solve for the vector difference between the velocities found in steps 1 and 2. 

(4) Find the velocity the vehicle has at the final point in the transfer ellipse. 



* = \Hr-^ orVl 



Vi Ti 

r« 
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(5) Compute the velocity required to keep the vehicle in the final orbit. 




(6) Find the vector difference between the velocities found in steps 4 and 5. 



(7) Find the total Av for the maneuver by adding the Av irom step 3 to the Av from 
step 6. 




Figure 37 



From a practical point, once the required Av is known, the amount of propellant 
required for the maneuver can be computed from 

v AY - = In = in mass ratio 

I.pg V W ) 

as discussed in Chapter 3. 

There are, of course, other ways to accomplish an altitude change. One such 
method is the Fast Transfer, useful when time is a factor. 

In the Fast Transfer, the transfer ellipse is not cotangential to the final orbit 
but crosses it at an angle (Fig. 37). Again the Av is applied in two increments, 
but Av 2 , applied at the intersection of the transfer ellipse and the target orbit, 
must achieve the desired final velocity in the proper direction . The steps for 
calculating the required Av are similar to those for the Hohmann transfer, noting 
that velocity differences must be treated as vectjrial quantities. For the same 
altitude change the fast transfer requires more Av. 

Other methods of changing altitudes are discussed in other texts on astro- 
nautics. The procedures are similar to those discussed here. The magnitude and 
direction of the velocity vector remains the critical factor. 



Plane Change 



Changing the orbital plane of a satellite also requires the expenditure of 
energy. This is apparent if the vector diagram representing two circular orbits 
of the same altitude is examined, the only difference being in their inclinations: 




Figure 38 

In this case, the orbital speeds, V] and v 2 , are equal except that they arc 90 
to each other. Transferring from the polar orbit to the equatorial orbit would 
require the Av represented oy the dashed arrow. For a 90° plane change (an 
extreme case) the Av exceeds the existing orbital speed. It should be noted that 
a change from one plane to another can only be accomplished at the intersection 
of the two planes. 

The Av required to change the orbital plane any spe fied amount can be deter- 
mined by examining the vectors involved. The problem is solvable by use of 
the Law of Cosines. For example, if v, represents the existing orbital velocity, 
v 2 the final orbital velocity, and a the desired plane change angle,* then: 




Av 



2 



2 2 o 

V -f V - ZV V COS oC 
1 2 12 



Figure 39 

If only the plane is to be changed, then v, = v 2 and the problem is simplified. 
But, if altitude (or eccentricity) is also to be changed, Vi and v 2 will not be equal. 

The amount of Av required to accomplish a plane change is, of course, de- 
pendent upon the amount of change desired. It is also a function of the altitude 
at which the change is made. 



• See appendix B. 
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Less Av is required to make a plane change at high altitudes than at low 
altitudes because the orbital speed of the vehicle is less at higher altitudes. In 
other words, it is more economical, in terms of propellant required, to make 
plane changes where the speed of the satellite is low — at apogee, or at high 
altitudes. 



Combined Maneuvers 

If a requirement exists to perform both a plane change and an altitude change, 
some economy will result if the operations are combined. 

The problem of combining a plane and altitude change is solved quite simply 
by considering the vector diagram. For example, it is desired to change the altitude 
of a vehicle from 100 NM circular orbit to 1500 NM circular orbit with a plane 
change of 10 3 . Recognizing that the plane change is more economically made at 
altitude, the plan is to combine the plane change with injection from the Hoh- 
mann transfer ellipse into the 1500 NM circular orbit. A typical Kohmann altitude 
change is initiated at a point of intersection of the two planes by increasing the 
vehicle’s speed. 

At the apogee of the transfer ellipse (1500 NM altitude) the vehicle’s speed 
is 19,800 ft/sec. The required circular speed is 21,650 ft/sec. The complete prob- 
lem, at apogee, looks like this: 

19,800 ft/ sec Speed at apogee of transfer ellipse 



Av Req’d 



Figure 40 




The Av required for the combined maneuver is calculated by use of the Law of 
Cosines and is Av = 4,055 ft/sec. It is also necessary to compute the angle at 
which the Av is to be applied. This calculation may be accomplished with the 
Law of Sines. 



PERTURBATIONS 

When the orbit of an artificial satellite is calculated by use of the assumptions 
given thus far, it will vary slightly from the actual orbit unless corrections are 
made to take care of outside forces. These outside forces, known as perturbations, 
cause deviations in the orbit from those predicted by two-body orbital mechanics. 
In order to get a better understanding of how the satellite’s actual orbit is going 
to behave, one must consider the following additional factors: 

( 1 ) The earth is not the only source of gravitational attraction on the satellite 
since there are other gravitational fields (principally of the sun and moon) in 
space. This effect is greatest at high altitudes (above 20,000 NM). 

2-50 




73 



(2) The earth is not a spherically homogeneous mass but has a bulge around 
the equatorial region. This additional mass causes the gravitational pull on the 
satellite not to be directed toward the center of the earth. This is the major 
perturbation effect at medium altitudes (between 300 NM and 20,000 NM). 

(3) The earth has an atmosphere which causes drag. This effect is most 
significant at low altitudes (below 300 NM). 



Third Body Effects 

One cause of perturbations is the introduction of one or more additional 
bodies to the sys. em creating a problem involving three or more bodies. ^gure 
41 shows an exaggerated perturbation, referred to a? a hyperbolic encounter. 
Initially, the satellite is in orbit 1 about an attracting body such as the earth. 




Figure 41. Hyperbolic encounter. 



As the satellite approaches the moon, the gravitational influence of the moon 
dominates, and the center of the moon becomes the focus instead of the center 
of the earth. Since the vehicle approaches the moon with more than escape 
velocity, it must leave the moon’s sphere of gravitational influence with greater 
than escape velocity. This means that the vehicle’s velocity with respect to the 
moon is greater than that required to escape. Therefore, for the s ort time a 
the moon is the attracting body, the vehicle is on a hyperbolic path with respect 
to the moon. When the satellite leaves the sphere of influence of the moon, 
switches back to the earth’s sphere of influence and goes into a new elliptical 
orbit about the earth. The next time the satellite returns to this region, the 
moon will have moved in its orbit, and the satellite, therefore, will now maintain 
orbit 2. A hyperbolic encounter is a method of changing the energy level of a 
satellite. By proper positioning, it could be used to increase or ecrease t e 
energy of a space vehicle. Of course, the change in energy of the vehicle is 
offset by the change in energy of the second body (the moon in the case 

illustrated) . 



Effects of Oblate Earth 



Another cause of perturbations is the bulge of the earth at the equator some- 
times called the earth’s oblateness. The effect of this oblateness on the satellite 
can be seen if we imagine the earth to be made up of a sphere which has an 
added belt of mass wrapped around the equatorial region. As shown in Fig. 42, 
the primary gravitational attraction F, directed to the center of the earth, will 
now be “disturbed” by the much smaller but nevertheless significant attractions 
Fi and F- directed toward the near and far sides of the equatorial bulge. With 
Tj smaller than r>, F, will be larger than F^, and the resultant force obtained 
by combining F, Fj, and Fu will now no longer point to the center of the earth 
but will be deflected slightly toward the equator on the near side. As the satellite 
moves in its orbit the amount of this deflection will change depending on the 
vehicle's relative position and proximity to the equatorial region. 




Figure 42. Deviation in force vector caused by the oblateness of the earth 

Two perturbations which result from this shift in the gravitational force are: 

( 1 ) Regression of the nodes. 

(2) Rotation of the line of apsides (major axis) or rotation of perigee. 

Regression of the nodes is illustrated in Fig. 43 as a rotation of the plane of 
the orbit in space. The resulting effect is that the nodes, both ascending and 
descending, move west or east along the equator with each succeeding pass. 
The direction of this movement will be opposite to the east or west component 
of the satellite’s motion. Satellites in the posigrade orbit (inclinations less than 
90°) illustrated in Fig. 43 have easterly components of velocity so that the 
nodal regression in this case is to the west. The movement of the nodes will be 
reversed for retrograde orbits since they always have a westerly component of 
velocity. Figure 44 shows why, for a vehicle traveling west to east, the regression 
of the nodes is toward the west. The original track from A to B would cross 
the equator at fi]. Simplifying the effect of the equatorial bulge to a single 



0 
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Figure 43. Regression of the nodes. 

impulse at point E, the track is moved so that it crosses the equator at fl 2 . At 
point F the simplified effect of the bulge is a single impulse down, changing 
the orbital path line to the line FD, which would have crossed the equator at ft 3 . 
This effec*, regression of the nodes, is more pronounced on low-altitude satellites 
than high-altitude satellites. In low altitude, low inclination orbits the regression 
rate may be as high as 9° per day. Fig. 45 shows how the regression rate 
changes for circular orbits at various altitudes and inclination angles. Note 
that nodal regression is zero in the polar orbit case. It has no meaning in 
equatorial orbits. 




Figure 44. Regression of nodes toward the west when vehicle is traveling west to east. 




Nodal regression rate per day for circular orbits. 
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Figure 45. 



Satellites requiring sun synchronous orbits (for photography or other reasons) 
are an example of how regression of the nodes can be used to practical ad- 
vantage in certain situations. In Fig. 46 the satellite is injected into an orbit 
passing over the equator on the sunlit side of the earth at local noon (the sun 
overhead). This condition initially aligns the orbital plane so that it contains a 
line between the earth and sun. The altitude of the near circular orbit deter- 
mines the angle of inclination required in order to maintain the sun synchronous 
nodal regression rate of approximately one degree per day (360°/year). In 
Fig. 45 note the required inclination angles of approximately 95° to 105° for the 
altitudes shown. If inclination angle and orbital altitude have been chosen 
correctly, then regression of the nodes will rotate the plane of the orbit (change 
the angle of right ascension) through 90° every three months as shown in Fig. 46. 
Thus, we see that this perturbing phenomenon due to the earth’s equatorial 
bulge maintains the desired angle of right ascension which, along with the angle 
of inclination, orients the orbital plane within the celestial sphere. Therefore, as 
the earth moves in its orbit the desired orientation for best photography is main- 
tained without using propellant. 




Figure 46. Sun-synchronous orbit. 
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Figure 47. Earth's equatorial bulge changes the argument of perigee. 

Rotation of the apsidal line is shown in Figure 7 (The apjdtd linens *e 
line joining apogee and perigee-the major axis I . : e cause of P" 
tion due to oblateness is difficult to visualize; however, ‘heresult is 
trajectory rotates within the orbital plane about the occupied focus. Th 
has the effect of shifting the location of perigee, thus changing the argument of 

Th" rTof "change in the argument of perigee is a function of satellite 
altitude and inclination angle. At inclinations of 63.4 and 116.6 the 
rotation is zero. Figure 48 illustrates how the apsidal rotation rate vanes with 
inclination angle for orbits with a 100 NM perigee and different apogee altitudes. 




Figure 48. Apsidal rotation rate per day far orbits with 100 NM perigee altitude. 
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Drag Effects 



Drag on a satellite will cause a decrease in eccentricity, a decrease in the 
major axis, and a rotation of the apsidal line. Figure 49 illustrates the change 
in eccentricity and major axis. The original orbit is represented by the curve 
Ai Ao Ai, with the focus at F, and a major axis equal to line segment Ai A 2 , 
Assuming that the drag is concentrated near perigee, the speed at Ai eventually 
will diminish to circular speed, and the new path will be Ai A s A u with the 
major axis decreasing to line segment At A 3 . After the orbit lias decayed to 
approximately zero eccentricity (a circle), further decay will result in a nearly 
circular spiral with ever-decreasing radius. 

Some other causes of perturbations are electromagnetic forces, radiation pres- 
sures, solar pressures, and gas-dynamic forces. 




Figure 49. Decrease in the eccentricity of a satellite orbit caused by drag. 

THE DEORBITING PROBLEM 

The general operation of moving a body from an earth orbit to a precise point on 
the surface of the earth is a difficult problem. TTu orb-/ may be circular or ellipti- 
cal, high or low, and inclined at various degree * * v rotational axis of the earth. 
The deorbiting maneuver may be the lofted, deprc.. c/, or retro approach which is 
the only method discussed in this section. Time oi flight is a major parameter in 
deorbiting and presents one of the more difficult of mathematical problems. But, the 
theory developed here is sufficiently general in scope that it may be modified to 
solve most problems. The complexity of the problem is reduced in the following sim- 
plified illustration of deor biting from a circular, polar orbit from a point over the 
North Pole, by assuming no atmosphere.* The approach and solution will permit 
impacting any earth target, provided fuel for retrothrust is no limitation. 



* The student may consider the effect of the atmosphere by using the radius of reentry rather than the 
radius of the earth, and by then considering the range and time of flight for a specific reentry body. The 
solution may also be modified to consider elliptical orbits and inclined orbits. 
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( '■sider a target at 60 north latitude. 11 the earth could be stopped fiom lotat- 
ing (see Fig. 50)*, with the target in the plane of the circular orbit, then it only 
would be necessary to apply a rctrovelocitv (Av,) to the orbital body so it would 
traverse the ellipse shown (are 1-4) and impact the target. The reentry velocity 
magnitude ( Vj> ) , equal to the circular speed (v,.) minus the reti ovelocity (Av,), 
would be necessary at apogee ( 1 ) to impact any target on the 60 north parallel. 
But, the direction of v t> would have to be oriented with respect to the earth’s axis so 
that it would lie in the plane formed by the earth's axis and the predicted position of 
the target. 

Since the earth is rotating at a constant speed, the time of flight (t b ) of the orbital 
object must be known in order to predict where the target will be at impact. For ex- 
ample, when the orbital vehicle is at Point 1 over the North Pole, the target is at 
Point 2. During the descent of the payload, however, the target moves to Point 3. In 
order to predict Point 3, t b must be known. 

Then, if the magnitudes of v e , v b , and a, the angle between these two vectors, are 
known, the Law of Cosines can be used to determine Avo. Application of Avo to v 0 
insures that v b will have the proper magnitude and direction so that the object will 
impact a selected target on a rotating earth. 



Looking first at the geometry of the problem (Fig. 51 ), notice that the total prob- 
lem lies in one plane. Basically, the radius of the circular orbit (r c ) and the latitude 
of the target (L) are known. Since the original circular orbit and the bombing trans- 
fer ellipse are coplanar and cotangential, the problem begins as a Hohmann transfer. 
To determine Av, the rearward velocity increment which must be applied to cause 
the object to impact the target canjpe computed from Av = v c — v b , if v c and v b 
are known. 



Write the general equation of the conic as k€= r (1 + € cos v ) and evaluate at 
the two known points, r,.** and r n , noting that at impact cos v = — cos 0t 
= — sin L, and at retrofire cos (180°) = — 1. 

Then ke = r 0 ( 1 — e b sin L) = r n ( 1 — e h ) 



Deorbiting Velocity 







* 8 figures related to the deorbit problem appear at the end of this section. 

* • r e is radius of the earth, but radius of re-entry could be used if desired. 
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Solving for e b . 



e b (r a ~ r 0 sin L) = r„ — r e 



, _ r a - r v 

€ b — ; f— 

r„ — r 0 sin L 

Note that in order to find e b and, in turn, the magnitude of v b , only the latitude 
of the target and the altitude of the circular orbit need be known. Suppose there is a 
satellite in a 500 NM circular polar oribit. When the orbital vehicle arrives directly 
over the North Pole, it is desired to deorbit an object which will impact at 60° north 
latitude. Assuming that the earth has no atmosphere and is nonrotating, what retro- 
vclocity is necessary? 



r„ - r p _ (23.94 X 10° ft) - (20.9 X 10° ft) 

r n - r e sin L (23.94 X 10 6 ft) - (20.9 X 10 6 ft) sin 60° 



c b = .520 


v l, — ^ -y— (1 — e b ) — -1 


/ 14.08 X 10 15 ft 3 /'sec 2 (1 - .520) 


V 23.94 X 10° ft 


v b = 16,790 ft/sec 





Av - v c — v b = 24,240 ft/sec — 16,790 ft/sec 
Av = 7,450 ft/sec 



Looking again at the geometry, note that there is a satellite in a 500 NM circular, 
polar orbit with v,. = 24,240 ft/sec. A letrovelocity increment, Av, equal to 7,450 
ft/sec was applied. This provided a magnitude v b = 16,790 ft/sec so that the object 
now follows arc 1-2 and impacts on 60° north latitude. 

Figure 52 provides a graph of velocity versus latitude to determine the magnitude 
of v b . The graph is for a specific altitude, for release from over the North Pole, and 
for no atmosphere. Included also are more general graphs, Figures 53 and 54. 



Deorbit Time of Flight 

As the bomb falls from apogee, the target moves toward the east due to earth ro- 
tation. Its speed is 1520 cos L ft/sec, so the necessity of accurately computing the 
time to bomb, t b , is readily apparent. 

First, look at the problem in schematic, Fig. 50. Recall that at the time the vehicle 
is over the North Pole at Point 1, the target is at Point 2. Thus, t b must be known so 
that the location of Point 3 can be predicted. If the meridian that Point 3 will be on, 
and the meridian with which the circular orbit coincides at the instant of deorbit are 
known, then the angle a can oe determined by subtraction. 

There are several methods for computing t b from release to impact. Eqn 10, App 
D could be used to determine t^, the time of flight for a ballistic missile launched 



from 6CT N latitude, with apogee at 500 NM over the North Pole, and impacting at 
60° N latitude. This would be the t^ to traverse arc 5-1-4 (Fig 50). The actual time 

to bomb from apogee would be half this amount. The general time of flight method 
is presented here. 



jj*- (u» - esinu,) - (uj - esinui) ( App D, Eqn 7 ) 



Since position 1 corresponds to apogee : 

Ui * 77, sin ni = 0 



Then t b = t,— 2 = \l~r ( u t “ «t» sin u t - tt) 

\ /X 

a = __Jjs 

1 4 €b 

_ e b — cos fl t — gb ~ sin L 
cos u t - j — e cos e t l - e 5 sin L 

All of these parameters are familiar except u which is the eccentric anomaly (Fig. 
55). If a perpendicular is dropped through the target to the major axis of the ellipse, 
it intersects a circle (with the center at C, and the diameter equal to the major axis 
of the ellipse) at Point Q. By definition, angle BCQ is u, the eccentric anomaly. The 
radius of the circular orbit, r a , is given. Both v b and e b can be calculated from for- 
mulas given previously. 

Working with the same example that was used to illustrate deorbiting velocity, t b 
will be calculated from only two known quantities — the altitude of the satellite and 
the latitude of the target. A satellite is in a 500 NM circular, polar orbit. Find the 
thne of flight, t b , from directly over the North Pole to a target on the 60° north par- 

al!- 1 



__ r„ — r e 


(23.94 X 10 6 ft) - (20.9 X 10 6 ft)_ 


€b r a - r e sin L 


“ (23.94 X 10° ft) - (20.9 X 10 6 ft) sin 60° 


e b = .520 




- - VC (1 - 


1 14.08 X 10 15 ft 3 /sec 2 (1 - .520) 
6 b) ~ \ 23.94 X 10® ft 


v b = 16,790 ft/sec 

r a _ 23.94 X 10® ft _ 15 7 o x 10 e f t 


* 1 + € b 
_ e b — sin L . 


1.520 

.520 - sin 60° _ _ ^ 


* 1 — e b sin L 


” 1 - .520 sin 60° 
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Noting u, lies between 180° and 270°: 



Ut = 360° - 129.6° = 230.4° = 4.084 radians 



tb 






— e b sin u t — 



■w] 



tb 



I (15.78 X 10° ft) 3 
\ 14.08 X 10 15 ft 3 /sec 2 



[4.084 - .520 sin 230.4° - 3.1416] 



t b = 676 sec = 11.27 min 



Figure 56 shows t b in minutes versus a plot of target latitude. Once again it must 
be recognized that this chart is for a specific orbital altitude and no atmosphere. See 
also the more general graph, Figure 57. 

Consider impacting a target at 60° north latitude and 30° east longitude. The 
plane of the satellite is coincident with the 70° east meridian. This means that, when 
the satellite is at Point 1 and the target is at Point 2, the angle between the target 
and orbital plane is 40°. However, the target is moving. It moves at: 



(360°/24 hrs) - ^ r - - = .25 deg/min 
60 min 

6 t = (11.27 min) (.25 deg/min) = 2.82° 
a — 40° — 2.82° = 37.18° 

Now apply the Law of Cosines and determine the Av 2 which must be applied to 
v c in order to impact the target: 



v_ = 24,240 a/sec 






Av 






/>N / N a= 37.13° 

v, = 16,790 f't/ficc 



VT= 






S = 43.15° 






/ ^ 




AVp = l4,8G0 ^t/sec 



Av 2 = (v b 2 + v c 2 — 2v b v c cos a) 1/2 

Av 2 = [(16,790 ft/sec) 2 + (24,240 ft/sec) 2 - 2(16,790 ft/sec X 24,240 ft/sec 

cos 37.18 °)] 1/2 

Av 2 = [2.82 X 10 8 + 5.88 X 10 8 - 6.49 X 10 s ] 1/2 
Av 2 = [2.21 x 10 8 ] 1 / 2 = 14,860 ft/sec 
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Also, the value of angle /? must be known so that the proper direction of Av 2 may 
be determined. The Law of Sines is preferred for this calculation, although the Law 

cf Cosines can be used. 



sin j3 
sin 



v b sin a 
Av 2 

16,790 ft/sec sin 37.18° 
14,860 ft/sec 



sin (3 = .6827 
/3 = 43.15° 



From the solution of this simplified object-from-orbit problem, it is apparent that 
such a calculation is not really simple. Other mathematical approaches and other op- 
erational problems are even more difficult. However, the theory presented can be 
extended to more difficult cases. 



Fuel Requirement 



It is also of interest to determine the propellant necessary to perform this maneu 
ver. Assuming an I sp = 450 sec (a reasonable figure in the near future) and an in 
itial weight of 10,000 pounds, compute the amount of propellant required, W p : 



In 



\ _ Av 14,860 ft/sec 

WT / “ U? (450 sec) (32.2 ft/sec 2 ) 



1.025 



= 10,000 jbs_ _ 3 590 lbs 

2.79 

= Wi - W 2 = 6,410 lbs 

This weight of propellant represents ^ = 64<1% of the weight m orblt 

prior to maneuvering. 

By restricting the plane change (side range) to a = 7.18 degrees, Av 2 = 7,875 

ft/sec, W p = 4,190 lbs, and = 41.9%. 

In summary, recall that the altitude of a satellite in a circular orbit and a time 
when the satellite was over the North Pole were given. A target was selected, and no 
earth atmosphere was assumed. From the theory, the required velocity for a deor lt- 
ing object was calculated, and the position of the target at the time of impact was 
predicted. With this information and the use of the laws of sines and cosines, the 
magnitude of retrovelocity and the direction to deorbit on target were calculated. 
Also, the amount of propellant required was computed. 
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2-62 




85 



1 




/ 







deorbit velocity 



I 



Figure 51 



2-63 

86 

O 

ERJC 



NORTH LATITUDE (L) IN DEGREES 




2-64 

87 

1 



ERIC 



NORTH PO 




sadjgdQ ui (*©) a-[Suv aSJu'eu 



2-66 



ERIC 



89 



Velocity increment (AV) in thousands feet/ second SOUTH POLE 







TIME OF FLIGHT GEOMETRY 



A 




Figure 55 



2 - 6 ? 




90 



NORTH LATITUDE (L) IN OEGREES 




TIME OF FLIGHT (t t ,) IN MINUTES 

Figure 56 




ERIC 



/I 






' vi 





S2 







